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Supplemental Declaration of Dr. James W. Mover under 37 C.F.R. § 1.132 

I, James W. Moyer, do hereby declare and state as follows: 

1. My credentials were presented in a previous declaration submitted May 
23, 2005. 

2. I am a named inventor on U.S. Patent Application No. 09/912,072 
(hereinafter "the '072 application"). 

3. At the time the work in the '072 application was carried out, AFLP 
technology had not been used to measure genetic diversity in poinsettia or in any 
species closely related to poinsettia. The application of AFLP technology to any 
particular plant species is considered uncertain and in light of this, to assess the 
relevance of the Loh et al., Barcaccia et al., Sukhwinder et al. and Barker et al. 
publications to poinsettia, it is further relevant to consider how distantly related 
poinsettia is to rice, willow, Pelargonium and Caladium. The taxonomic relationship of 
poinsettia as compared with these plants is shown in Appendix A. 

4. All of the plants at issue are in the Division Magnoliophyta, the division to 
which all flowering plants belong. Two of the plants referenced, rice and caladium, are 
further classified as monocots (Class Liliopsida), while three, pointsettia, willow and 
geranium, are dicots (Class Magnoliopsida). 

5. Rice is further classified into the Subclass Commelinidae, Family 
Poaceae, and Genus Oryza. The Commelinidae contains an enormous diversity of 
plants including rushes and sedges, cereals, pineapples, ginger, and arrowroot. 
Therefore, it is nearly impossible to generalize findings with regard to one member of 
the Commelinidae such as a grass (e.g., maize or rice) to another member such as a 
bromeliad (e.g., a pineapple). 

6. Caladiums are in the Subclass Aracidae. The Aracidae includes the 
palms, duckweed and caladium. Caladium is further classified into the Family Araceae 
(aroid), and Genus Caladium. The aroid group includes edible plants as well as those 
that are poisonous and includes philodendron, calla lily, Mexican breadfruit and taro. 
Thus, there is an enormous diversity of plants in the Subclass Aracidae and it would be 
difficult to generalize regarding one member such as duckweed with another member 
such as palm or philodendron. 

7. Within the class Magnoliopsida, or the dicots, willow (Salix) is further 
classified into the Subclass Dilleniidae. The Dilleniidae is once again a very diverse 
grouping including, for example, the mustard family (e.g., broccoli), the heath family 
(e.g., rhododendron and blueberry) and the cucurbit family (e.g., cucumber, squash). 
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Willow is further classified into the Family Salicaceae and the genus Sa//x. Once again, 
it would be very difficult to generalize findings about a member of this diverse Subclass, 
such as willow or broccoli, to another member such as blueberry or squash. 

8. Both Pelargonium and poinsettia are in the Subclass Rosidae. This is an 
incredibly diverse group of plants with 108 families and 58,000 species. Some 
examples of plant families that fall within the Rosidae include the carrot family, the 
apple family, the legume family (e.g., pea), and the dogwood family. In addition, the 
Rosidae subclass also includes the geranium family to which Pelargonium belongs and 
the very large and diverse euphorbia or spurge family (approximately 300 genera and 
7,500 species) to which poinsettia belongs. (Park et al., Int J. Plant Sci. 161:425-434 
(2000); C.L. Porter, Taxonomy of Flowering Plants, W.H. Freeman & Co., 472 pp., 
p.338, (1967)). Pelargonium is further classified into the Genus and Species 
Pelargonium peltatum. Poinsettia is in the Genus Euphorbia, which has been 
described as one of the largest and most complex genera of flowering plants with about 
1600 species. Id.. Clearly, it is nearly impossible to generalize findings about one 
member in the Subclass Rosidae, such as carrot to another member such as dogwood. 
This would be difficult even within the single Genus Euphorbia, of which poinsettia is a 
member, due to its great size and diversity. 

9. Plant scientists interested in determining genetic relationships between 
and within plant cultivars, varieties or species would be well-aware of the distant 
relationship between poinsettia and the plants studied in the cited publications, and 
would not find our achievement in poinsettia obvious in light of results in a tree such as 
willow, a cereal plant such as rice, an aroid plant such as caladium, or even geranium. 
It is simply not the case that AFLP analysis in any other plant, particularly distantly 
related plants such as those in the cited publications, would suggest application of 
AFLP analysis in poinsettia. A scientist in this field would at most think to try such a 
technique but would enter into the research without any expectation of success. This 
would be particularly true in the case of poinsettia since it is known to have a very 
narrow genetic base (See, prior Moyer Declaration submitted May 23, 2005). 

10. I hereby declare that all statements made herein of my own knowledge 
are true and that all statements made on information and belief are believed to be true; 
and further that these statements were made with the knowledge that willful false 
statements and the like so made are punishable by fine or imprisonment, or both, under 
Section 1001 of Title 18 of the United States Code and that such willful false statements 
may jeopardize the validity of the application or any patent issued thereon. 
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APPENDIX A. Taxonomic relationships between poinsettia and the plant species in the cited 
references. 



Kingdom: Planta 

Subkingdom: Tracheobionta - vascular plants 
Superdivision: Spermatophyta - seed plants 

Division: Magnoliophyta - flowering plants 

Class: Magnoliopsida - dicotyledons 

Subclass: Dilleniidae 
Order: Salicales 

Family: Salicaceae 

Genus: Salix - Willow 

Subclass: Rosidae 

Order: Geraniales 

Family: Geraniaceae 

Genus: Pelargonium - Geranium 

Order: Euphorbiales 

Family: Euphobiaceae 

Genus: Euphorbia - Poinsettia 

Class: Liliopsida - monocotyledons 

Subclass: Arecidae 
Order: Arales 

Family: Araceae 
Genus: Caladium 

Subclass: Commelinidae 
Order: Cyperales 
Family: Poaceae 

Genus: Oryza- Rice 
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A PHYLOGENETIC STUDY OF TRIBE EUPHORBIEAE (EUPHORBIACEAE) 

Ki-Ryong Park 1 and Wayne J. Elisens 

/ Department of Biology, Kyung-Nam University, Masan 631-701, Korea; and Department of 

Botany and Microbiology, University of Oklahoma, Norman, Oklahoma 73019, U.S.A. 

A phylogenetic investigation of a monophyletic lineage of spurge plants, tribe Euphorbieae, was conducted 
to elucidate evolutionary relationships, to clarify biogeographic patterns, and to reexamine the previous clas- 
sification of Euphorbieae. Cladistic analyses of the 52 morphological characters of 61 species resulted in 2922 
equally most parsimonious trees of 193 steps with a consistency index of 0.34. The strict consensus tree 
indicates genus Anthostema of subtribe Anthosteminae as a likely sister group to all other members of tribe 
Euphorbieae. The morphological data support a monophyletic origin of subtribe Euphorbiinae, but the sub- 
tribes Anthosteminae and Neoquillauminiinae did not form monophyletic groups. Although the previous 
taxonomic treatments within tribe Euphorbieae have supported the generic status of Pedilanthus, Monadenium y 
Synadenium, Cha?naesyce, and Elaeophorbia, the results of this analysis do not support generic placement of 
them based on cladistic principles. Recognition of these groups as genera results in Euphorbia becoming a 
paraphyletic group. One solution to this problem in Euphorbieae is to divide the largest genus Euphorbia 
into several monophyletic genera and to keep the generic ranks for previously recognized genera. The distri- 
bution of basal endemic genera in Euphorbieae showed African and east Gondwanan affinities and strongly 
indicated that the ancestor of Euphorbieae originated prior to the breakup of Gondwanaland from an old 
group in Euphorbiaceae. However, some recent African taxa of Euphorbia should be interpreted by transoceanic 
dispersal from the New World ancestors. 
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Introduction 

The reconstruction of phylogenetic relationships of organ- 
isms facilitates the testing of systematic hypotheses as well as 
organizing of species into a formal classification (Nelson and 
Platnick 1981; Funk and Brooks 1990; de Queiroz and Gau- 
thier 1994). The phylogenetic relationships of tribe Euphor- 
bieae have yet to be thoroughly examined using cladistic meth- 
ods (Gilbert 1994; Webster 19946). The diversity of the 
taxonomic treatments within this tribe is cumulatively a result 
of the tribe's large size, the profusion of intergrading and over- 
lapping characters (Sherff 1940), and a tribal classification that 
occurred without a comprehensive study of the entire group 
and throughout the entire range of its species. Thus, testing 
of generic and infrageneric boundaries within tribe Euphor- 
bieae based on cladistic principles is one of the fundamental 
tasks of concern to systematists who specialize in Euphorbia 
and Euphorbiaceae. 

Tribe Euphorbieae (ca. 2000 species), one of the most dis- 
tinctive and largest tribes in the Euphorbiaceae (Croizat 1938), 
is characterized by the possession of a unique cyathium that 
consists of a central carpellate flower and four or five groups 
of basal male-flower clusters (Croizat 1937; Dressier 1957; 
Webster 19946). The tribe is widely distributed throughout 
the world, but the majority of species occurs in tropical and 
subtropical Africa and America. The most recently proposed 

' Author for correspondence; e-mail park@kyungnam.ac.kr. 
Manuscript received March 1999; revised manuscript received October 1999. 



classification of tribe Euphorbieae (Webster 1975, 19946) rec- 
ognizes three subtribes: Anthosteminae, Neoguillauminiinae, 
and a true cyathium-bearing Euphorbiinae. This classification 
has been based mainly on staminate calyx, petaloid appendage, 
and gland characters. Pax and Hoffman (1931) recognized 10 
genera and represented their relationships with a phyletic 
graph (reconstructed in fig. 1; Pax 1924). In this phyletic graph, 
Pax recognized four groups and proposed that Anthostema 
and Dichostemma are placed in a close relationship; Calyco- 
peplus was treated as the most closely related genus to a group 
of Euphorbia + Elaeophorbia, and Synadenium, Stenadenium, 
and Monadenium share a common ancestor. Pax may have 
overemphasized the shape of involucral glands, and the genera 
with continuous glands (i.e., Synadenium^ Monadenium, and 
Stenadenium) were treated as an independent lineage — one dif- 
ferent from the true cyathium-bearing Euphorbia in his 
diagram. 

An alternative phyletic scheme, one that includes Neoguil- 
lauminia and Pedilanthus, was proposed by Croizat (1937), 
and it has contributed significantly to our understanding of 
the phylogenetic relationships and cyathium evolution of tribe 
Euphorbieae (fig. 1). Croizat (1937) noted the presence of a 
male or female calyx in Anthostema, Dichostemma , Neoguil- 
lauminia, and Calycopeplus and segregated them from the true 
cyathium genera that have glands attached to the rim of the 
involucre. He considered Calycopeplus to be closely related to 
Euphorbia. Additionally, he emphasized the zygomorphic cy- 
athium of Cubanthus and Pedilanthus and treated them as 
groups outside of Euphorbia. 
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Fig. 1 Relationships within Euphorbieae proposed previously by 
Pax (1924) and Croizat (1937). Hypotheses of relationships are por- 
trayed cladistically. 

Within Euphorbia s.l., the New World endemics Chamae- 
syce and Pedilanthus have been treated as distinct genera or 
as closely related subgenera of subg. Agaloma based on their 
shared petaloid appendages, which appear around the invo- 
lucral glands (Wheeler 1941; Dressier 1961). The species of 
Chamaesyce differ from the rest of the Euphorbiinae in that 
they have several unique characters, and they have been rec- 
ognized as a genus by some authors (Croizat 1943; Webster 
19946). The polyphyly of Euphorbia s.l. was questioned based 
on the diverse origin of glands in Euphorbieae (Croizat 1938). 

Cytological evidence (Hans 1973) indicates that within Eu- 
phorbia s.l., modern euphorbia species evolved from the prim- 
itive dendroid subg. Esula (x=6, 7, 8, 9, 10) (Perry 1943) or 
subg. Euphorbia (x=10) through an aneuploid or polyploid 
series (Webster 1967; Mehra and Choda 1978). Anthostema 
is distinguished from the rest of Euphorbieae by a distinct base 
number x=ll. Common basic chromosome number (n=14) 
indicates a close relationship among native New World subg. 
Agaloma, Poinsettia, and Chemaesyce (Urbatsch et al. 1975). 

A recent cladistic analysis of 33 morphological characters 
at the sectional level was conducted for New World subtribe 
Euphorbiinae (Park 1996). However, this phylogenetic analysis 
of New World Euphorbiinae did not support the generic rec- 
ognition of Chamaesyce because it was nested within subg. 
Agaloma and because it required generic segregation for many 
groups within Agaloma (Park 1996). This analysis also did 
not strongly support the monophyletic origin of New World 
endemic petaloid-appendaged taxa and nested a monophyletic 
Chamaesyce, Poinsettia, and Pedilanthus within a paraphyletic 
Agaloma. 

Using morphological characters, we conducted a phyloge- 
netic investigation of spurge plants, tribe Euphorbieae, to elu- 
cidate phylogenetic relationships among the major taxonomic 
groups of Euphorbieae and to clarify biogeographic patterns 
among them. Additional objectives were to test the monophyly 
of subgenera proposed by Wheeler (1943) and by Carter and 
Radcliffe-Smith (1988) in Euphorbia s.l. and to examine 
whether New World petaloid-appendaged taxa form a mono- 
phyletic group. 

Material and Methods 

Monophyly, Terminal Taxa, and Outgroup Selection 

The Euphorbieae have been accepted as a monophyletic 
group that is characterized by its unique cyathium (Croizat 



1937; Webster 1967; Mehra and Choda 1978; Gilbert 1994). 
Recent phylogenetic analysis of subfamily Euphorbioideae 
strongly supported a monophyletic tribe Euphorbieae (K.-R. 
Park and A. Backlund, unpublished data) and its sister group 
relationships to species in tribe Hippomaninae + Hureae. 
Based on these analyses, two species of Hippomaninae, Ma- 
prounea brasiliensis and Sebastiana brasiliensis, were selected 
as outgroups for the ingroup taxa. 

Although it is ideal to start a cladistic analysis with mono- 
phyletic species groups as terminal taxa in the analysis of large 
ingroups (Sanderson 1991), this was not possible in the case 
of Euphorbieae because no prior cladistic hypotheses of re- 
lationships within Euphorbieae were present, except for New 
World taxa. Thus, species were selected as terminal taxa, and 
at least one species of most higher taxonomic groups in Eu- 
phorbieae was chosen (Boissier 1862; Pax and Hoffman 1931; 
Carter and Radcliff-Smith 1988; Webster 19946) in this anal- 
ysis (appendix). Fifty-nine ingroup species were selected as 
terminal taxa. 

Characters and Their States for Cladistic Analysis 

Fifty-two qualitative characters are used in this analysis (ta- 
ble 1). Characters were obtained by direct examination of her- 
barium specimens and from the taxonomic literature. 

Habit. In the Euphorbieae, most species outside of Eu- 
phorbia possess a woody habit with many stems (character 1). 
Within Euphorbia s.l., however, various forms of habit are 
present, even within the same subgenus. Trees or shrubs are 
quite common in subg. Euphorbia, whereas the woody con- 
dition is usually restricted to species occurring in insular or 
coastal habits in subg. Chamaesyce and Esula (Park 1996). 
Succulent stems (character 2) are very common in African en- 
demic members of subg. Euphorbia, Lacanthis, Eremophyton, 
and genus Elaeophorbia; succulence is generally associated 
with deciduous leaves (character 12). With ternate stems and 
opposite phyllotaxy, stem articulation (character 3) occurs only 
in Agaloma sect. Alectoroctonum (Park 1996). 

Inflorescence branching. The floral rays (character 7) of 
herbaceous species of Euphorbiniiae usually branch dichasially 
several times, and they are subtended by an equal number of 
leafy bracts. Some subtropical members of Euphorbiniiae pos- 
sess a pseudoumbel dichasium at the tips of main stems. An 
umbel-dichasial inflorescence is present in most Euphorbia spe- 
cies; a thyrse (character 8) occurs only in some species in subg. 
Agaloma sect. Alectoroctonum and in Dichostemma. The pres- 
ence of scattered rays (character 9) is considered to be ancestral 
within Euphorbieae, based on the previous analysis (Park 
1996), but condensed or simple rays are restricted to Eu- 
phorbia subg. Poinsettia, a few species of subg. Agaloma, subg. 
Tirucalli and Trichadenia. 

Leaves, bracts, and spines. The leaves (character 11) of 
Euphorbieae (in general) are alternate, whereas those of species 
of sect. Zygophyllidium in Agaloma and Poinsettia are op- 
posite. Species of Agaloma sect. Alectoroctonum (Euphorbia 
cotinifolia and Euphorbia leucocephala) have verticillate 
leaves. Leaf margins of most species of subg. Eremophyton, 
Poinsettia, and Chamaesyce have serrate or dentate leaves 
(character 13), whereas remaining species in Euphorbieae have 
entire leaves as an ancestral condition. The equilateral bases 
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Table 1 



Characters and Character States Usee 


in Cladistic Analysis of Tribe Euphorbieae (Euphorbiaceae) 


Character 




State 


1 Wokit- 


0 


= tree or shrub; 1 = perennial herb; 2 = annual herb 






= woody; 1 = semisucculent; 2 = succulent 




0 


= nonarticulate; 1 = articulate 




o 


= absent; 1 = present 




A 
u 


= absent; 1 = present 




n 
u 


= cylindrical; 1 = angled 






= absent; 1 = present 




o 


= absent; 1 = present 




o 


= scattered; 1 = condensed; 2 = simple 




o 


= absent; 1 = present 




o 

u 


= alternate; 1 = opposite 




A 

u 


= persistent; 1 - deciduous 




n 

u 


= entire; 1 = dentate or serrate 




o 


= equilateral; 1 = nonequilateral 




o 


= absent; 1 = present 




o 


= absent; 1 = present 






= absent; 1 = present 




A 
u 


= inclined (elliptic); 1 = even (half-circle shape or deltoid) 






= green only; 1 = green and anthocyanic 




A 


= absent; 1 - present 




U 


= absent; 1 = present 




r\ 
U 


= terminal; 1 = axillary 




n 


= absent; 1 = present 




u 


= enclosed by bract; 1 = not enclosed by bract 


-) r \ X „ 1 ~ (\ 


n 
u 


= unorganized; 1 = group of four; 2 = group of five 


26. Staminate calyx 


u 


= present; 1 = absent 


27. Carpellate flower 


A 

u 


= basal; 1 ~ apical 




A 

u 


= present; 1 = absent 




A 

u 


= smooth; 1 = verrucose 




A 

u 


= absent; 1 = present 




A 

u 


= three-lobed or rounded; 1 = triangular 




n 

U 


= straight; 1 = reflexed 




A 

u 


= undivided; 1 = divided 




A 
\J 


= connate all of the length; 1 = free; 2 = connated at base 




A 

u 


= absent; 1 = present 






= round; 1 = deltoid; 2 = divided at margin; 3 = pectinate with acute tips 


■J —i ^l__J „L„„- 


A 


= rounded; 1 - crescent 




A 
U 


= four; 1 = five; 2 = one or two; 3 = more than five 




A 
U 


= free; 1 = connate with ring shape 


40 Gland surface 


0 


= flattened or folded; 2 = funnel- or cup-shaped 


41. Gland position 


0 


= between staminate flower clusters; 1 = between ihvolucral lobes 




0 


= absent; 1 = 4-Iobed; 2 = unlobed 




0 


= not petaloid; 1 = petaloid 




0 


= globose; 1 = angular 


45. Seed dorsal view 


0 


= rounded; 1 = conical; 2 = rectangular 




0 


= present with round shape; 1 = present with cap shape; 2 = absent 


47. Seed tubercula 


0 


= absent; 1 = present 


48. Seed honeycomb 


0 


= absent; 1 = present 


49. Seed horizontal ridges 


0 


= absent; 1 = present 




0 


= thick; 1 = thin 




0 


= round; 1 = deeply three-lobed 




0 


= more than one; 1 = one 



of leaves are plesiomorphic in Euphorbieae, whereas the non- 
equilateral leaves (character 14) with interpetiolar stipules 
(character 15) are restricted to species of Euphorbia subg. 
Chamaesyce (Croizat 1943; Webster 1994b). Only Monaden- 
ium, some species of subg. Esula, and Tirucalli have semicir- 



cular or deltoid bracts, but in all other taxa, elliptic bases of 
bracts are predominant (character 18). Species having succu- 
lent stems with spiny outgrowths (character 21) are restricted 
to the African endemic subg. Euphorbia, subg. Lacanthis, and 
Elaeophorbia, The spine shield (character 20) surrounding the 
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base of leaves is an advanced character of subg. Euphorbia in 
Euphorbia s.l. (Carter 1994). 

Staminate and carpellate flowers. The staminate calyx 
(character 26) of Anthostema, Diehostemma, and outgroup 
taxa of Euphorbieae is present, but it is lacking in most of the 
taxa within Euphorbieae. Subsequently, Webster (19946) rec- 
ognized subtribe Anthosteminae, including the above two gen- 
era, based mainly on this character. In some species of Cu- 
banthus and Euphorbia subg. Agaloma, the vestigial forms of 
the staminate calyx apparently are present (Croizat 1937; Park 

1998) , but in this analysis, they are coded as lacking the sta- 
minate calyx. A pistillate calyx (character 28) is found in sub- 
tribes Anthosteminae, Neoguillauminiinae, and only Cuban- 
thus of Euphorbiinae. All other taxa have no pistillate sepals, 
but they are well developed in outgroup taxa. 

Involucral characters. Most of the members of the New 
World subtribe Euphorbiinae are characterized by petaloid ap- 
pendages on the involucral glands (character 35). Wheeler 
(1939), Dressier (1957, 1961), and Webster (1967) hypothe- 
sized close relationships among Euphorbia subg. Chamaesyce, 
subg. Agaloma, and Pedilanthus based on this character. Two 
different gland positions (character 41) are known in Eu- 
phorbieae. The glands of subtribe Euphorbiinae are on the 
margin of the cyathia, which is opposite to the clusters of male 
flowers, whereas the glands of the two remaining subtribes are 
inserted between the clusters of male flowers. The glands of 
outgroups, subtribe Hippomaninae, occur on each side of the 
floral bracts enclosing the male flowers. Therefore, the al- 
ternative position of the glands in Anthosteminae and Neo- 
guillauminiinae is a symplesiomorphic condition within the 
Euphorbieae. 

Cladistic Analysis 

Multistate characters were treated as unordered (Fitch 
1971). The data sets were analyzed using maximum-parsimony 
algorithms available on the PAUP version 4.0b2 (Swofford 

1999) for the Macintosh. The matrix was analyzed using 
random taxon addition sequences followed by the efficient 
tree bisection-reconnection branch-swapping algorithm with 
MULPARS. A strict consensus tree was generated. In order to 
evaluate the stability of different branches in the obtained trees, 
a bootstrap analysis with^ 1000 replicates of data (using fast- 
heuristic search, MacClade version 3.0; Maddison and Mad- 
dison 1992) was used to explore character evolution and bio- 
geographic pattern and to depict alternate tree topologies. An 
area cladogram was produced by replacing the name of taxa 
on the data matrix with their geographic occurrences (Platnick 
and Nelson 1978). 

Results 

Cladistic analyses resulted in 2922 equally most parsimo- 
nious trees of 193 steps with a consistency index of 0.34 and 
a retention index of 0.74. The strict consensus tree is repre- 
sented in figure 2. 

The strict consensus tree indicates genus Anthostema of sub- 
tribe Anthosteminae as a likely sister taxon to all other mem- 
bers of tribe Euphorbieae. The monophyly of subtribe Eu- 
phorbiinae, a true cyathium-bearing group, was supported by 



three synapomorphic characters, but subtribes Anthosteminae 
and Neoguillauminiinae did not form a monophyletic group 
(fig. 2). All trees supported sister-group relationships between 
subtribe Euphorbiinae and Calycopeplus oligandrus of Neo- 
guillauminiinae and between Cubanthus linearifolius and the 
remaining taxa of Euphorbiinae. Within subtribe Euphorbi- 
inae, the segregate genera Pedilanthus, Monadenium, Syna- 
denium, and Elaeophorbia were nested within a paraphyletic 
Euphobia s.l. 

The Monadenium + Synadenium clade forms a sister-group 
relationship with the Australian endemic Euphorbia subg. Er- 
emophyton clade, including Euphorbia parvicaruncula, E. 
stevenii, E. tannensis, and E. boophthona (fig. 2). Subgenus 
Eremophyton is paraphyletic. 

The strict consensus tree indicated a sister-group rela- 
tionship between Elaeophorbia and Euphorbia subg. La- 
canthis, which in turn are sister to Euphorbia subg. Eu- 
phorbia. The New World endemic Pedilanthus is sister to 
the clade consisting of Elaeophorbia, Euphorbia subg. La- 
canthis, and subg. Euphorbia. 

This analysis supports the monophyly of five subgenera 
within Euphorbia s.l.: Chamaesyce, Poinsettia, Tirucalli, Eu- 
phorbia, and Lyciopsis (fig. 2). Three subg. Agaloma, Ere- 
mophyton, and Esula are found to be polyphyletic or para- 
phyletic. Only one additional step was required for subg. Esula 
to become monophyletic. 

Discussion 

Phylogeny and classification. A recent classification sys- 
tem (Webster 19946) in tribe Euphorbieae recognized three 
subtribes and initially segregated Anthosteminae from Neo- 
guillauminiinae + Euphorbiinae based on the presence of a 
staminate calyx in the Anthosteminae subtribe. Croizat (1937) 
suggested that Anthostema had basal placement within the 
tribe and implied that it was not closely related to Dieho- 
stemma. This analysis supported Croizat's hypothesis and in- 
dicated that the staminate calyx is the plesiomorphic condition 
within the tribe. In our analyses, monophyly of the Anthos- 
teminae (Pax 1924; Webster 1994b), Anthostema and Dieho- 
stemma, required only one additional step, and this subtribe 
cannot be fully rejected based on this analysis. 

Croizat (1937) also hypothesized a close relationship be- 
tween Euphorbia and Calycopeplus and proposed Calycope- 
plus as an ancestor of Euphorbia s.l. This analysis supported 
Croizat's (1937) hypothesis, if Cubanthus is considered a 
member of Euphorbia s.l. Croizat treated Cubanthus and Pedi- 
lanthus as an independent lineage within the Euphorbieae, but 
our data and analysis did not support this relationship. In this 
study, Pedilanthus was nested within Euphorbia s.l. and 
formed a sister group to the African endemic clade consisting 
of Elaeophorbia, Euphorbia subg. Lacanthis, and Euphorbia. 
This result represented a departure from previous hypotheses 
related to the origin of Pedilanthus. Dressier (1957) and Park 
(1996) hypothesized that Pedilanthus was derived from the 
New World endemic subg. Agaloma sect. Tricherostigma. 
However, early interpretations of relationships (Boissier 1862; 
Millspaugh 1913; Croizat 1942) indicated a close affinity be- 
tween Cubanthus and Pedilanthus. 

Among genera of Euphorbieae, Pax (1924) proposed a close 
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.■■ Anthosteminae 
mm Neoguillauminiinae 
w Euphorbiinae 



S. brasiliensis, M. brasiliensis 

A. aubrganum I 

D. gfaucescens I 
N. cleopatra I 
C oligandrus 1 
C linearifolius £ 
£ drupifera $ 
£ gottlbie $ 
£ eyassiana $ 
£ robecchii % 
P. tithymaloides $ 
£ marginata $ 
M. iaeve £ 
S. compactum $ 
£ parvlcaruncu!a$ 
£ stevenii $ 
£ tannensis $ 

E. boophthona $ 
£ polyantha $ 
E. etuberculosa $ 
£ longituber. £ 
£ systyia £ 
£ tophiosperma £ 
E. agowensis $ 
£ pirottae % 
£ chersonea % 
£ heterophyfla $ 
£ dentata % 
£ sarcodes % 
£ hexagonoides £ 
£ exstipulata jj 
E. hexagona $ 
£ bilobata $ 
E. uniglandulosa $ 
£ francoana $ 
£ multiseta $ 
£ ocyrnoidea £ 
£ graminea £ 
£ atoto £ 
£ supina £ 
£ angusta £ 
£ ariensis $ 
£ innocua $ 
E. corollata Z 

y 

£ cuneata £ 

£ yoyae ^ 

£ scheffleri % 

£ grantii % 

E. cafifornica $ 

£ misera $ 

E. leucocephata $ 

£ cotini folia $ 

£ nubica $ 

£ tirucalli £ 

£ sieboldiana % 

£ pekinens'ts $ 

£ ebracteolata % 

£ dendroides £ 
£ plumeroides 



Outgroup 
>4nt/josre/na 
Dichostemma 
Neguiflauminia 
Calycopeplus 
Cubanthus 
Eiaeophorbia 
subg. Lacs n this 
subg. Euphorbia 
subg. Euphorbia 
Pedilanthus 
subg. Agaloma 
Monadenium 
Synadenium 
subg. Eremophyton 
subg. Eremophyton 
subg. Eremophyton 
subg. Eremophyton 
subg. Eremophyton 
subg. Eremophyton 
subg. Eremophyton 
subg. Eremophyton 
subg. Eremophyton 
su bg . Eremophyton 
subg. Eremophyton 
subg. Agaloma 
subg. Poinsettia 
subg. Poinsettia 
subg. Agaloma 
subg. Agaloma 
subg. Agaloma 
subg. Agaloma 
subg. Agaloma 
subg. Agaloma 
subg. Agaloma 
subg. Agaloma 
subg. Agaloma 
subg. Agaloma 
subg. Chamaesyce 
subg. Chamaesyce 
subg. Chamaesyce 
subg. Agaloma 
subg. Agaloma 
subg. Agaloma 
Lyclopsls 
Lyciopsis 
subg. Lyciopsis 
subg. Trichadenia 
subg. Agaloma 
subg. Agaloma 
subg. Agaloma 
Agaloma 
Tirucalli 
subg. Tirucalli 
subg. Fsu/a 
subg. Esu/a 
subg. Esu/a 
subg. Esu/a 
subg. Esu/a 
subg. Esu/a 



subg. 
subg. 



subg. 
subg. 



Fig. 2 Strict consensus tree of 2922 equally most parsimonious cladograms of tribe Euphorbieae. Genera and subgeneric names of Euphorbia 
are represented. Numbers below branches denote bootstrap values above 50%. 



relationship between Synadenium and Monadenium and seg- 
regated them accordingly from Calycopeplus, Euphorbia, and 
Eiaeophorbia. Our cladistic analysis supports this relationship 
and indicates common ancestry between these genera, but our 
analysis did not support segregation of the genera from £«- 
phorbia, because they showed a close affinity with the Ew- 
phorbia subg. Eremophyton group. 



Several workers recognized Eiaeophorbia Stapf as a genus 
separate from Euphorbia, based on its unique indehiscent dru- 
paceous fruit (Brown 1915; Pax and Hoffmann 1931; Hutch- 
inson 1969; Carter and Radcliffe-Smith 1988). Recently, Evans 
and Kinghorn (1977), using diterpene profiles, rejected the 
recognition of Eiaeophorbia. Our analysis of morphological 
characters placed Eiaeophorbia as a sister group to Euphorbia 



430 



INTERNATIONAL JOURNAL OF PLANT SCIENCES 




A. aubrganum 

D. glaucescens 
N. deopatra 
C oligandrus 
C linearifolius 
£ drupifera 
£ gotttbie 

E. eyassiana 
£ robecchii 
P. tithymaloides 
E. marginata 
M. taeve 
S. compactum 
£ parvicaruncula 
£ stevenii 
E. tannensis 
£ boophthona 
E. polyantha 
E. etubercufosa 
£ longituber. 
E. systyla 
E. lophiosperma 
E. agowensis 
E. pirottae 
E. chersonea 
£ heterophylla 
E. dentata 
£ sarcodes 
E. hexagonoides 
E. exstiputata 
E. hexagona 
E. uniglandulosa 
E. bilobata 

> E. francoana 
» E. multiseta 

> E. ocymoidea 

* E. graminea 

* £ atoto 

> £ supina 
1 £ angusta 
i £ ariensis 

> £ innocua 
i £ corollata 
. £ caiifonica 

> £ misera 
£ ieucocephata 

> £ cotinifolia 
£ cuneata 
£ joyae 

> E. scheffleri 
t £ grantii 

£ nubica 
£ tirucalli 
£ sieboldiana 
£ pekinensis 
£ jolkini 
£ ebracteolata 
£ dendroides 
£ piumeroides 



Fig. 3 One of the most parsimonious cladograms from morphological analysis and the geographic distribution of each species {1 - America; 
2 = Africa; 3 = Australasia; 4 = Asia). Closed circles represent species that have appendaged glands on the rims of the involucre. 



subg. Lacanthis, one that was nested within the subg. Eu- 
phorbia + Lacanthis clade. 

Within the genus Euphorbia, subg. Euphorbia is one of the 
most well recognized groups with a unique spine shield and 
pair of spines (Carter 1994). Gilbert (1987) and Carter (1994) 
hypothesized that subg. Euphorbia may have originated in 
Asia from a subg. Tirucalli-Wke ancestor. Carter (1994) also 
suggested that the calloused leaf scars of subg. Tirucalli should 



be homologous with the spine shield in subg. Euphorbia. The 
results of this study do not support the relationships proposed 
by Gilbert (1987) and Carter (1994), but they do strongly 
indicate the close relationships among subg. Euphorbia, La- 
canthis, and Elaeophorbia. 

Euphorbia subg. Eremophyton (Boiss.) Wheeler is a group- 
ing of ca. 20 species native to Africa and Australia, and it is 
usually divided into two or three sections based on seed char- 
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acters (Boissier 1862; Hassall 1977; Carter and Radcliffe- 
Smith 1988). Carter and Radcliffe-Smith (1988) recognized 
only two sections: sect. Pseudacalypa, with a conical, ecarun- 
culate seed, and sect. Eremophyton, with a cap-like carun- 
culate seed. The presence of several independent lineages 
within subg. Eremophyton had been proposed by several au- 
thors (Boissier 1862; Croizat 1972; Radcliffe-Smith 1974). 
This study did not support the monophyly of subg. Eremo- 
phyton. However, the Australian endemic species of subg. Er- 
emophyton formed a monophyletic group, which is in accord 
with the hypothesis proposed by Hassall (1977). 

A close relationship among New World petaloid-append- 
aged groups of Euphorbiinae has been proposed by several 
monographers (Wheeler 1939; Dressier 1957, 1961; Webster 
1967), but close relatives of the New World taxa were not 
identified. Our analysis indicated that the taxa with petal oid 
appendages were paraphyletic and shared a common ancestry 
with African endemic subg. Trichadenia and subg. Lyciopsis 
(fig. 3). If this relationship is supported, then the finger-like 
processes on the glands in the latter two subgenera are ho- 
mologous to the petaloid appendages on the glands in New 
Worid taxa. 

The close relationships of Euphorbia subg. Chamaesyce and 
the probable derivation from New World taxa have been in- 
dicated by the distribution of several characters: abortion of 
the main stem, a C 4 photosynthetic pathway, interpetiolar stip- 
ules, and nonequilateral leaves (Koutnik 1987). Wheeler 
(1941), Hurusawa (1954), and Webster (1967) hypothesized 
several models for the origin of Chamaesyce from Euphorbia 
subg. Agaloma, whereas Croizat (1943) insisted on an inde- 
pendent origin of Chamaesyce from Euphorbia. A recent clad- 
istic analysis of New World Euphorbiinae (Park 1996) sup- 
ported Norton (1900) and Webster's (1967) views that 
Chamaesyce is closely related to subg. Agaloma species such 
as Euphorbia ipechacuanhae and Euphorbia innocua. The re- 
sults of our analysis strongly supported the origin of Chamae- 
syce within subg. Agaloma. 

Dressier (1961) postulated that Euphorbia subg. Poinsettia 
was closely allied to subg. Agaloma sect. Dichillum Boiss. Park 
(1996) hypothesized that the condensed inflorescence of sect. 
Petaloma might be homologous with that of subg. Poinsettia, 
whereas Croizat (1942) speculated on a close relationships 
with African and Australian subg. Eremophyton species such 
as Euphorbia crotonoides and Euphorbia eremophila. This 
analysis, which includes Old World Euphorbieae, did not sup- 
port a close relationship between subgenera Poinsettia and 
Agaloma but indicated subg. Poinsettia nested within the Eu- 
phorbia subg. Eremophyton clade, close to Euphorbia agow- 
ensis and Euphorbia pirottae. 



Within a polyphyletic Euphorbia subg. Agaloma, this study 
also supported the monophyly and sister-group status of sec- 
tions Alectroctonum (Euphorbia cotim folia and Euphorbia 
leucocephala) and Tricherostigma (Euphorbia calif ornica and 
Euphorbia miser a). 

Although previous taxonomic treatments within tribe Eu- 
phorbieae have recognized the generic status of Pedilanthus, 
Monadenium, Synadenium, Chamaesyce, and Elaeophorbia, 
the results of this analysis do not support generic recognition 
of these taxa. Recognition of these groups as genera results in 
Euphorbia becoming a paraphyletic group. 

Biogeography. The Euphorbieae, the largest tribe in Eu- 
phorbiaceae, is broadly divided into two geographical areas: 
a Paleotropical distribution in primitive taxa, and a Neotrop- 
ical and Eurasian distribution in derived members. Most au- 
thors strongly advocate a recent origin for Euphorbieae and 
accept long-distance dispersal to explain the worldwide pat- 
terns within the tribe (Webster 1994a). In this study, the dis- 
tribution of basal endemic genera in Euphorbieae showed Af- 
rican and east Gondwanan affinities (fig. 3) and strongly 
indicated that the ancestor of Euphorbieae originated prior to 
the breakup of Gondwanaland from an old group in 
Euphorbiaceae. 

The pattern within Euphorbia s.l. indicates that the origin 
and basal diversification occurred in Australasia and Asia and 
was followed by dispersal to Africa and the New World. Eu- 
phorbia subg. Tirucalli consists of ca. 30 species and shows 
disjunct distributions through the Paleotropical region in Mad- 
agascar, the Cape region, East Africa, and Indochina. The clad- 
istic trees revealed that the species of the subg. Tirucalli, one 
of the primitive members of Euphorbia s.l., originated from 
the common ancestor with subg. Esula in Southeastern Asia 
and then dispersed to India and moved on to diversification 
in Africa (ca. 20 species) from Tertiary. The most recent African 
members of Elaeophorbia and Euphorbia (subg. Eremophy- 
ton, subg. Euphorbia, subg. Lacanthis) are independently in- 
troduced from the New World. Although the pattern observed 
in basal groups in Euphorbieae may be explained by Gand- 
wanan vicariance, some recent Euphorbia s.l. should be in- 
terpreted based on transoceanic dispersal models (Webster 
1994a). 
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Appendix 



Table A1 

Data Matrix of Character States for 52 Morphological Characters Used in 
Cladistic Analysis of Euphorbieae 

1111111111222222222233333333334444 444 444555 
12345678 9012345678 9012345 678 90123456789012345678 9012 

Tribe Euphorbieae: 

Subtribe Anthosteminae: 
Genus Anthostema (2): 

A. aubrganum 0000000010000000000001101010000010070100010000000000 

Genus Dicbostemma (3): 

D. glaucescens 00000001 00000000000000101 010010010070000010007000001 

Subtribe Neoguillauminiinae: 

Genus Neoguillauminia (1): 

N. cleopatra 0000000000000000000000101110000010070300021000000001 

Genus Calycopeplus (5): 

C. oligandrus 0000010010110000000001101110000010070000121000000101 

Subtribe Euphorbiinae: 
Genus Cubanthus (3): 

C. linearifolius 00000010100000000000001121100001 12070202120072000101 

Genus Elaeophorbia (1): 

E. drupifera '.. 0200010000010000000011112111000010070100120002000001 

Genus Monadenium (50): 

M. laeve 1 10000000000000 10100011 121110001 10070210120121 100101 

Genus Synadenium (10): 

S. compactum 0100000000000001000000112111000110070210120007100101 

Genus Pedilanthus (15): 

P. titbymaloides 0100000010010001000001112111000112100000120002000101 

Genus Euphorbia (1000): 

E. dendroides 0000001000000000110000112111000110071000120000000101 

£. plumerioides 00000010000000001000001121110001100 70100120000000101 

E. sieboldiana 11000010000000001100001121110001100 71000120000000101 

E. pekinensis 1100001000000000100000112111100110070000120000000101 

E. jolkini 1100001000000000100000112111100110070000120000000101 

£. ebracteolata 1100001000000000110000112111000110070000120000000101 

£. eyassiana 0201010010010001000111112111000010070100120002100171 

£. robecchii 0201010010010001000111112111000010070100120002000171 

E. gottlbie 0201110000010001001011112111000710070100120002000171 

£. nubica 0001001-020010001110000112111000110070000120000000171 

E. tirucalli 00010010200100011100001121110001100 70100120000000101 

E. cuneata 000110100 7000001100000112111010010070101120002000171 

E.joyae 0000001007000001100000112111010010070100120002000171 

E.scheffleri 0000001020000001100000112111000010130100120002000171 

E. grantii 0000001000000001110000112111000110130000120002000171 

E. etuberculosa ..... 1200001000000001100000112111000110070201120112001101 

£. longituberculosa 1200001007000001100000112111000110070201120112001101 

E. systyla 2100001007000001100000112111000110070000120112101101 

£. polyantha 1100001000000001110000112111000110070000120121100101 

£. lophiosperma 1100001000001001100000112111010110070000120112101101 

£. agowensis 2100001001001001100000112111010111070000120121100101 

E. pirottae 2100001001000007100000112111010111070000120121100101 

£. parvicaruncula 2100000001011001000000112111000110070000120121100101 

E. stevenii 2100000001011001000000112 1110001100 70000120002 1 00101 

£. tannensis 2100000001011001000000112111000110070000120121100101 

£. boophthona 2100000001011001000000112111000110070000120121100101 

E. innocua 1100007000000001700000112111010110100000120102010101 

E. corollata 1100001000000001100000112111000110100100120002010101 

E. hexagonoides 2100000007100001000000112111000111110100120112100171 

E. bilobata 2100000007100001000000112111010111120100120112010171 

E. exstipulata 2100000007101001000000112111010111100000120110100171 

E. hexagona 2100000007101001000000112111000111110100120000100101 

£. californica 0001100020000000000000112111000110100100120002010171 
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Table A1 



(Continued) 

111111111122222222223333333333444 4 44 44 4 4555 
1234 567 89012345678901234 56789012345678 9012345678 9012 

E. misera 0001100020000000000000112111010110100100120002010171 

E. leucocephala 0010000100200001101000112111000110100100120100010171 

E. cotinifolia 00100001002000001010001121110001101001001201020101?1 

£. chersonea 2100001010001001100000112111010111 ??00011201121001?1 

E. uniglandulosas ...... 21000000010000010000001121110?01111201001201121101?1 

E. graminea 11000000070000010000001121110001 11 1000001201121101 ?1 

E. ariensis 1100001000100001101000112111000111100000120112010171 

E. francoana 2100000001000007000000112 1110101 1 11201001201120101 ?1 

£. multiseta 21000000010000010000001121110?01101201001201120101?1 

£. ocymoidea 2100000001000000000000112111010111120000120112010101 

£. marginata 2100001010000001101000112111010110100100120002100101 

£. sarcodes 0000001000000007100000112111000110120000120112100101 

£. heterophylla 2100000010101000101000112111000111070201120112100101 

£. dentata 2100000010101007101000112111000111070201120110100101 

£. atoto 1100000001?00111000000112111000111100000120002000111 

£. supina 2100000001?01111000000112111011111100000120112001111 

£. augusta 1100000001700111000000112111011111100000120112001111 

Tribe Hippomaninae: 

M. brasiliensis 0000000070000000000000000000000100070700000000000000 

S. brasiliensis 0000000070000000000000000000000000070700000000000000 



Note. "?" designates missing states; classification system based on Webster (1994fc). Total number of 
species within the genus is given in parentheses. 
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Part III I SELECTED ORDERS AND FAMILIES OF DICX>TYLEDOl|l 



ORDER EUPHORBIALES 

An anomalous group, including only the following family. Hutchinson^ 
suggests a possible origin from Malvales and Sapindales. The plants seenjjL 
to be highly specialized, as evidenced by very reduced flowers in many j^ 
the peculiar and flower-like inflorescence in Euphorbia, and the milky juice |jj 
and attendant complex secretory tissue. 

EUPHORBIACEAE: Spurge Family 
Herbs, shrubs, or trees, often with milky juice. Leaves simple or compound, | 
usually alternate, sometimes reduced to spines. Flowers monoecious or || 
sometimes dioecious, commonly in cymes, with or without a perianth, 
sometimes with a corolla. Staminate flowers variable, often reduced to a .|g| 
single stamen. Pistillate flowers rather uniform, consisting of a single pistil ^ 
of mostly 3 carpels, 3-lobed, and forming a capsule or schizocarp which splits |p 
into three 1-seeded nutlets in fruit. In specialized forms, such as Euphorbia, 
the inflorescence may simulate a flower, being reduced to one or more || 
cymules or cyathia, these sometimes with ornamental bracts. ||f 
The family includes about 283 genera and 7,300 species, widely distributed , || 
in temperate and tropical regions, the Indo-Malayan region and Brazil 
being the chief centers of distribution. Some African species of Euphorbia .;|| 
resemble cactus plants. Economically important products include food, 
drugs, and rubber. Many species are poisonous, the genus Toxicodendron of 
South Africa including some of the most poisonous plants known. 

examples: Euphorbia: the largest genus, with 1,600 species estimated, 
including E. pulcherrima, the Poinsettia. 
Hevea braziliensis: an important rubber tree. 
Manihot utilissima: Cassava, furnishing tapioca and arrowroot 

starch from the tuberous roots. 
Ricinus communis: Castor Bean, the source of castor oil. 



-Wr^ /axcrooTTY^ Vjauxn t^s t HTSLp^, vm*. Freeman 
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Abstract Amplified Fragment Length Polymor- 
phism (AFLP) analysis was used to generate informa- 
tive DNA fingerprints from Magnolia denudata, M. 
liliiflora, and putative cultivars of their hybrid cross, 
i.e., M. x soulangeana, M. xsoulangeana 4 Lennei\ 
M. x soulangeana 'Rustica Rubra 1 , M. x soulan- 
geana 'Alba', M. x soulangeana Tluby\ and M. x 
soulangeana 'San Jose'. A single specimen of un- 
known origin and representatives of the closely re- 
lated Af. campbellii and M. stellata were also included 
in the analysis. We show how this method may be 
used for distinguishing among cultivars. 

Keywords Amplified Fragment Length Polymor- 
phism (AFLP); cultivar; DNA fingerprinting; 
hybrid; Magnolia 



INTRODUCTION 

A cultivar must be clearly distinct, uniform and sta- 
ble in its characteristics (Trehane et al. 1995). Dis- 
tinguishing among cultivars based on morphology 
can be difficult because the diagnostic characters are 
often small or may not be apparent all year round. 
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DNA fingerprinting provides a means to test the 
hybrid origin of cultivars and cultivar similarity and 
may lead to a better understanding of our cultivated 
plants (Lee etal. 1996). 

Magnolia x soulangeana or the saucer magnolia, 
resulted from a cross between A/, denudata and M. 
liliiflora made by Etienne Soulange-Bodin in 1820. 
Soulange-Bodin established and was the first direc- 
tor of the Royal Institute of Horticulture near Paris 
where he produced this hybrid (Spongberg in Hunt 
1998). Although commonly listed as M. x 
soulangiana, the correct spelling for this hybrid is 
with the termination -eana as described by Hunt 
(1998). Callaway (1994), considered it likely that 
Soulange-Bodin' s cross was also made independ- 
ently in Japanese nurseries, possibly occurring spon- 
taneously when the two parent plants were grown in 
close proximity. After the introduction of the hybrid 
by Soulange-Bodin, a plethora of named forms were 
introduced, including seedlings from backcrosses 
and open pollinations. These forms cover the com- 
plete range between the two parents in flower colour, 
shape, and size. It is now practically impossible to 
keep track of all of the forms and account for their 
origins. 

A rootstock of unknown origin in the Magnolia 
collection at Lincoln University, New Zealand that 
grew after the death of the scion M. acuminata 
'Golden Glow' produced flowers similar to M. x 
soulangeana cultivars in shape and colour, but were 
generally much larger. Flowering period and mor- 
phological characters were insufficient to distinguish 
this taxon, which has potential as a new cultivar for 
amenity horticulture. 

The DNA fingerprinting technique known as Am- 
plified Fragment Length Polymorphism (AFLP) 
(Vos et al. 1995; Mueller & Wolfenbarger 1999) 
may be applied directly to distinguish among and 
identify cultivars (Baraccia et al. 1998; Dirlewanger 
et al. 1998; Nadarajan et al. 1999; O'Hanlon et al. 
1999; Tignon et al. 2000). The use of DNA profil- 
ing for plant variety registration is attracting particu- 
lar attention (Law et al. 1998). The AFLP method 
allows the generation of 1 0- 1 00 times more markers 
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per reaction than other fingerprinting techniques, the 
large number of putative loci sampled providing a 
measure of variation across a wide portion of the 
genome (Sharma et al. 1996). 

In the present investigation, we used fluorescent 
AFLP analysis to study genetic variation in Af. 
denudata Desr., M. liliiflora Desr., and cultivars of 
their hybrid cross. Af. campbellii Hook. & Thorns, 
and Af. stellata (Sieb. & Zucc.) Maxim, were in- 
cluded for reference because these species are esti- 
mated (Qui et al. 1995) to be among the closest 
relatives of Af. denudata and M. liliiflora. The main 
question we addressed was: can Magnolia cultivars 
be distinguished using AFLP markers? 



MATERIALS AND METHODS 

Plant materials and DNA samples 

Samples from 12 mature Magnolia specimens were 
used for the present study (Table 1 ). A single cultivar 
of unknown origin was included for investigation, 
plus samples of M. campbellii and Af. stellata. Two 
leaves were used for separate extractions from M. x 
soulangeana and carried through the complete AFLP 
process as a control. All selective amplifications 
were repeated for the first primer combination and 
subsequently carried through as a test of reproduc- 
ibility. 

DNA extraction, digestion, and ligation 

Total DNA was extracted using a modification of the 
CTAB method (Doyle & Doyle 1987). DNA was 
precipitated in ethanol and resuspended in 25 [ll TE 



(10 mAf Tris-HCl, pH 8.0; 0.1 mAf EDTA) buffer. 
DNA concentrations were estimated and standard- 
ised against known concentrations of Lambda DNA 
on 1% agarose gels. 

Initially 10 ^1 of the total DNA extract was di- 
gested with 6 nl digestion buffer (50 mAf Tris-HCl, 
pH 7.5; 50 mAf (CH 3 COO) 2 Mg.4H 2 0; 250 mAf 
CH 3 COOK), 5 nl dH 2 0, 1 \il of each restriction 
enzyme; EcoRl (10 U/^l, Roche Molecular Bio- 
chemicals) and Tru9 1 (an isoschizomer for MseT) (10 
U/^l, Roche Molecular Biochemicals), for 2 h at 
37°C followed by 1 5 min at 70°C. A 5 |xl aliquot was 
run on a 1% agarose gel to check digestion was 
complete and the remaining DNA was then ligated 
by adding 1 \il of each of the EcoRl (5 pmol/|ll) 
and Msel (50 pmol/|al) adapters (Life Technologies), 
2 til of ligation buffer (100 mAf Tris-HCl, pH 8.3; 
15 mAf MgCl 2 ; 500 mAfKCl), 1 |xl T4 DNA ligase 
(1 U/jll, Gibco), and 5 ^1 dH 2 0. The mixture of 
ligation cocktail plus digested DNA was incubated 
overnight at 4°C. 

Pre-and selective amplifications 

The total pre-amplification reaction volume was 40 
This reaction consisted of 2 jxl ligated DNA (di- 
luted 1 :5), 4 |il of 10 x PCR reaction buffer (100 mAf 
Tris-HCl, pH 8.3; 15 mAf MgCl 2 ; 500 mAf KC1), 
0.8 fil 10 mAf dNTP, 31.04 nl dH 2 0, 0.16 nl Taq 
Polymerase (5 U/|il, Roche Molecular Biochem- 
icals), and 1 \il of each pre-amplification primer, 
each with one selective nucleotide; EcoRl+ 1 primer 
5'-G ACTGCGTACCAATTCA-3' mdMseH 1 5- 
G ATG AGTCCTG AGTAC-3' (50 ng/^1, Life Tech- 
nologies). Pre-amplification was carried out in a 



Table I Magnolia species, cultivars, and hybrids, Lincoln University, New 
Zealand, living plant and herbarium (LINC) voucher numbers for samples used 
in this study. 



Magnolia species/cultivar 

M. campbellii 

M. stellata 

M. denudata 

M. liliiflora 

M. liliiflora 'Nigra* 

M. x soulangeana 

M. x soulangeana 'Lennei' 

M. x soulangeana 'Rustica Rubra' 

M. x soulangeana *Alba' 

M. x soulangeana 'Ruby' 

M. x soulangeana 'San Jose* 

M. x "unknown" 



Plant no. 


Herbarium no. 


None 


LINC20004 


59 


LINC20007 


39 


LINC20010a 


62 


LINC20008 


3 


LINC20005 


27 


LINC20009 


46 


LINC20032 


1 


LINC20003 


14 


LINC20035 


78 


LINC20006 


2 


LINC20036 


18 


LINC20002 



I 
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thermocycler using the following temperature pro- 
file: 40 cycles of 30 s at 94°C, 30 s at 50°C, 60 s at 
72°C. 

The total selective amplification reaction volume 
was 20 |il and contained 1 \il of preamplification 
product (diluted 1 : 1 0), 2 ^1 1 0 x PCR reaction buffer 
without MgClj (as supplied with Expand High Fi- 
delity PCR system, Roche Molecular Biochemicals), 
2.5 nl 25 mA/MgCl 2 , 0.5 ^1 10 mA/dNTP, 1 1.8 pi 
dH 2 0, 0.2 |xl Taq Polymererase and 1 \il of each of 
the selective primers (10 ng/|il EcoRI and 50 ng/jil, 
Msel, Life Technologies). £coRI primers were la- 
belled with HEX. Selective amplifications were 
made using four primer combinations: (1) £coRl- 
AGC/A&el-CAG; (2) EcoRl-AGC/Msel-CAT, (3) 
£coRl-AGC/MseI-CTG; (4) EcoRX-ACGIMsel- 
CAT. The following cycle profile was used for se- 
lective amplification: 5 cycles of 30 s at 94°C, 30 s 
at 65°C, 60 s at 72°C, followed by 6 cycles of 30 s 
at 94°C, 30 s at 60°C, 60 s at 72°C, followed by 24 
cycles of 30 s at 94°C, 30 s at 56°C, 60 s at 72°C. 

Generation of AFLP fingerprints 

From the 20 |xl selective amplification product, 4 \il 
was run out on a 3% agarose gel to check for dis- 
tinct amplification products. The remaining 16 ^1 
was purified via ethanol precipitation and aliquots 
of the products run out on an ABI 377 automatic 
sequencer. Fragments were detected by laser and 
accurately sized with an internal standard. Digitally 
converted raw data were saved as samples migrated 
past the fluorescence detector. Multilocus profiles 
were visualised using ABI GENESCAN software. 
Fragment differences between samples were identi- 
fied as those differing by 1 bp. AFLP profiles were 
scored for presence/absence of fragments between 
50 and 400 bp and a binary matrix prepared. 

Data analyses 

The statistical software package SYSTAT vers. 5 
was used to calculate the percentage congruence 
between replicate samples using the simple match- 
ing coefficient. 

SYSTAT vers. 5 was also used to compare three 
levels of genetic variation: (1) between replicate 
selective amplifications; (2) between M. x soulan- 
geana cultivars; and (3) between species. A similar- 
ity matrix was generated for primer combination 1 
based on the coefficient of Jaccard (S) (Sneath & 
Sokal 1973). This coefficient does not consider 
the shared absence of a character as indicating 
similarity: 



S=a/(a + f> + c) 

where a = number of fragments present for both 
samples; b = number of fragments present for 
Sample A, but not for Sample B; c = number of 
fragments present for Sample B, but not for Sample 
A. When comparing variation between samples, 
the similarity between every possible pair was 
used to calculate the median similarity index (Sokal 
& Rohlf 1981) for each level of variation present 
among the operational taxonomic units (OTUs) in 
primer combination 1. 

Using the numerical taxonomy program NTS YS- 
pc (Rohlf 1990), similarity matrices were generated 
from each of the binary matrices for each primer 
combination based on Jaccard's coefficient. A 
pairwise correlation was calculated between the el- 
ements of each of these similarity matrices and the 
Mantel test (Mantel 1967), with 100 random permu- 
tations performed to test for significant similarities 
between primer combinations. 

Binary matrices for each primer combination 
were modified so that each contained only those 
OTUs shared among all matrices. The matrices with 
identical OTUs were then pooled to form a single 
matrix. The coefficient of Jaccard was then used to 
generate a matrix of similarities for. ( 1 ) cluster analy- 
sis using the unweighted pair-group method arith- 
metic average (UPGMA); and (2) principal 
coordinates analysis (PCOORDA). A UPGMA 
dendrogram was constructed and the first two co- 
ordinates were plotted from the PCOORDA. 

The program AFLPapp (Benham, J. 1997: 
AFLPapp. A tool for analysis of AFLP data: http:// 
hordeum.msu.montana.edu/software/software.html) 
was used to identify the fragments inherited from 
each putative parent for all primer combinations. The 
method takes the first parent and compares every 
fragment to the second parent. Fragments were not 
used in the analysis if: (1) a fragment from the first 
parent is also present in the second parent; and (2) 
if a fragment was not shared between any of the 
cultivars and parents. The Magnolia of unknown 
origin was included for analysis. 



RESULTS 

The congruence between replicate amplifications 
was 91% (81-100%). Median levels of similarity 
among OTUs in primer combination 1 based on the 
Jaccard coefficient were as follows: (1) between 
replicate selective amplifications, 85%; (2) between 
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Fig. 1 Plot of principal coordi- 
nate 2 versus 1 for Magnolia 
amplified fragment length poly- 
morphisms from four selective 
primer combinations. M.xsoulan- 
geana (a) and (b) represent differ- 
ent leaves from the same tree. 



Table 2 Polymorphism detected with four Amplified Fragment Length 
Polymorphism (AFLP) primer combinations for 13 Magnolia samples. 



AFLP primer 


No. of 


No. of polymorphic 


% 


combination 


fragments 


fragments 


polymorphism 


£coR1-AGC/Mk?I-CAG 


283 


225 


79.5 


EcoR\-AGC/Mse\-C AT 


447 


328 


73.4 


fcoRl-AGC/Afcei-CTG 


328 


276 


84.1 


£a>Rl-ACG/AfeeI-CAT 


343 


309 


90 


Total 


1401 


1138 


81.2 



Af. x soulangeana cultivars, 64%; and (3) between 
species, 55%. 

A total of 1401 AFLP fragments were recorded 
for 13 OTUs for the four pooled AFLP primer 
combinations, with an average polymorphism of 
8 1 .2% (Table 2). The Mantel test showed a highly 
significant (P - <0.001) correlation between all 
combinations of primer matrices. The matrix cor- 
relation statistics for each primer combination were: 
/• = 0.62 for primer 1 versus 2, r ~ 0.69 for primer 
I versus 3, /■ = 0.69 for primer 1 versus 4, r - 0.85 
for primer 2 versus 3, r = 0.86 for primer 2 versus 
4, r = 0.86 for primer 3 versus 4. Good separation 
was achieved using principal coordinates analysis, 
with the first coordinate explaining 64.2% of the 
variation. The second and third coordinates ex- 
plained only 5. 1 and 4.3% of the total variation re- 
spectively. Principal coordinate axes I and 2 (Fig. 
1 ) and UPGMA clustering (Fig. 2) were generally 



in concordance, both showing that: (1) M. 
campbellii and Af. stellata were separated from 
other OTUs; (2) M. x soulangeana cultivars 
grouped more closely with Af. denudata than with 
M. lililflora\ (3) samples originating from separate 
leaves of Af. x soulangeana grouped closely to- 
gether; and (4) the sample of unknown hybrid ori- 
gin appeared among the OTUs representing Af. 
denudata and A/, x soulangeana cultivars. No ties 
were encountered during UPGMA analysis. 

Magnolia x soulangeana cultivars contained 
greater proportions of fragments inherited from Af. 
denudata compared with M. liliiflora (Fig. 3). Num- 
bers of fragments shared between cultivars and par- 
ents ranged from Af. x soulangeana 'Alba' (165 
fragments shared with A/, denudata and 73 frag- 
ments shared with A/, liliiflora), to Af. x soulangeana 
'Lennei* (139 fragments shared with Af. denudata 
and 99 fragments shared with Af. liliiflora). 
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Fig. 2 Unweighted pair-group 
method arithmetic average 
(UPGMA) analysis of Magnolia 
amplified fragment length 
polymorphisms from four selective 
primer combinations. M. x 
soulangeana (a) and (b) represent 
different leaves from the same tree*. 
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Fig. 3 Number of amplified frag- 
ments from four primer combina- 
tions shared between Magnolia x 
soulangeana cultivars and their 
parents, M denudata (□) and M. 
Uliiflora (■). 
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DISCUSSION 

The results show AFLP data are useful for detect- 
ing variation within and among Magnolia cultivars. 
The 9 1% congruence for replicate AFLP profiles are 
comparable with previous studies (Huys 1996; Win- 
field et al. 1998) that reported 95-100% variation for 
reference samples, indicating the error present is 
similar to that commonly encountered The presence 
of error in the form of non-homologous AFLP prod- 
ucts, may bias genetic distance estimates between 
^axa, Rieseberg ( 1 996) suggested that Random Am- 
plified Polymorphic DNA (RAPD) homology was 
a function of taxonomic distance, i.e., the more 
closely the compared taxa are, the greater the prob- 
ability that a shared co-migrating band is homolo- 
gous. This might also be extended to the analysis of 



AFLPs. In the present study, the median similarities 
based on the Jaccard coefficient were generated 
specifically for the purpose of comparing different 
taxonomic levels. These showed decreasing percent- 
age similarity from replicate selective amplifications, 
to cultivars of M. x soulangeana and finally samples 
representing different species. 

Magnolia x soulangeana cultivars grouped most 
closely with M. denudata^ based on PCOORDA 
and UPGMA clustering (Fig. 1 and 2). All cultivars 
inherited a higher number of fragments from 
M. denudata compared with A/. Uliiflora (Fig. 3). 
We were able to estimate that the cultivar of un- 
known origin included in our study may have origi- 
nated from a cross between Af. denudata and M. 
Uliiflora, because it was most closely related to the 
M. x soulangeana cultivars 'San Jose 1 and 'Lennei' 
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(Fig. 1 and 2) and compared with the other cultivars, 
contained a similar number of fragments inherited 
fromM denudata and M. liliiflora (Fig. 3). Unpub- 
lished flowering records of magnolias at Lincoln for 
3 years from 1988 showed that the unknown mag- 
nolia was more similar to M. x soulangeana 'San 
Jose* than M. x soulangeana 'Lennei' in terms of the 
spring flowering period and accumulated thermal 
time. But M. "unknown" was unlike both M. x 
soulangeana 'San Jose' and M. x soulangeana 
'Lennei' in that only these two cultivars of M. x 
soulangeana have a second and significant flower- 
ing period in late summer. M x soulangeana 'Ruby' 
has also shown some tendency to flower in late sum- 
mer, but to a lesser extent. A significant period of 
flowering in late summer was also recorded for Af. 
liliiflora and M. liliiflora 'Nigra' for each of the 3 
years. 

Possible explanations for the close relationship 
implied for M. x soulangeana cultivars with M. 
denudata include: (1) the specimens do not repre- 
sent true M. x soulangeana, and are in fact simple 
cultivars of M. denudata; (2) the specimens are M. 
x soulangeana, but segregation has occurred in the 
F 2 and succeeding generations; and (3) the speci- 
mens are M. x soulangeana, but backcrossing to Af. 
denudata has occurred. Spongberg (1998) consid- 
ered many forms of M x soulangeana to approach 
or merge with one or other parent, possibly as a re- 
sult of segregation and introgression. Future inves- 
tigations of M. x soulangeana hybrid origin should 
include multiple representatives of each parent, par- 
ticularly M. denudata. 

AFLP variation may have potential for profiling 
Magnolia cultivars, providing evidence in addition 
to morphology, for registration purposes. However, 
it may be advantageous to apply selection criteria to 
fragments included for analysis, e.g., Escaravage et 
al. (1998), who considered only fragments of certain 
intensity in their analysis of clonal diversity in a 
Rhododendron population. In this way it may be 
possible to reduce error that might result from scor- 
ing artefact fragments that may result from partial 
digestion of the template genomic DNA, poor am- 
plification of fragments during PCR, etc. Tignon et 
al. (2000) reported uniformity of AFLP profiles 
among representatives of most apple cultivars they 
tested, but one cultivar showed variability between 
specimens, which indicated to them, a level of 
genetic instability within this cultivar. This high- 
lights the need for adequate testing of multiple 
representatives of each cultivar (and also from 
generation to generation) before the use of AFLP 



profiling as a tool for registration purposes. Accu- 
rate identification of cultivars would also require a 
database containing DNA fingerprints for all recog- 
nised cultivars for comparison. 
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Abstract DNA-based fingerprinting technologies have 
proven useful in genetic similarity studies. RFLP is still 
most commonly used in the estimation of genetic 
diversity in plant species, but the recently devel- 
oped PCR-based marker techniques, RAPDs, SSRs 
and AFLPs, are playing an increasingly important 
role in these investigations. Using a set of 33 maize 
inbred lines we report on a comparison of techniques to 
evaluate their informativeness and applicability for the 
study of genetic diversity. The four assays differed in 
the amount of polymorphism detected. The informa- 
tion content, measured by the expected heterozygosity 
and the average number of alleles, was higher for SSRs, 
while the lowest level of polymorphism was obtained 
with AFLPs. However, AFLPs were the most efficient 
marker system because of their capacity to reveal sev- 
eral bands in a single amplification. In fact, the assay 
efficiency index was more than ten-fold higher for 
AFLPs compared to the other methods. Except for 
RAPDs, the genetic similarity trees were highly corre- 
lated. SSR and AFLP technologies can replace RFLP 
marker in genetic similarity studies because of their 



Communicated by F. Salamini 

I. Pejic * V. Kozumplick 
Faculty of Agriculture University of Zagreb, 
Deptartment of Plant Breeding, Genetics and Biometrics, 
Svetosimunska 25, HR-10000 Zagreb, Croatia 

M. Morgante 

Dipartimento di Produzione Vegetale e Tecnologje Agrarie, 
Universita di Udine, Via delle Scienze 208, 33100 Udine, Italy 

P. Ajmone-Marsan • P. Castiglioni • M. Motto (S) 
Istituto Sperimentale per la Cerealicoltura, Via Stezzano 
24, 24126 Bergamo, Italy 

G. Taramino 

DuPont Agricultural Biotechnology, 

Delaware Technology Park, Suite 200, 1 Innovation Way, 

PO Box 6104, Newark, DE 19714-6104, USA 



comparable accuracy in genotyping inbred lines se- 
lected by pedigree. Bootstrap analysis revealed that, in 
the set of lines analysed, the number of markers used 
was sufficient for a reliable estimation of genetic sim- 
ilarity and for a meaningful comparison of marker 
technologies. 

Key words Zea mays L. * Genetic relationship • 
Molecular markers • DNA-fingerprinting • 
Genetic diversity 



Introduction 

Knowledge of germ plasm diversity and of relationships 
among elite breeding materials has a significant impact 
on the improvement of crop plants (Hallauer et al. 
1988). In maize, this information is useful in planning 
crosses for hybrid and line development, in assigning 
lines to heterotic groups, and in plant variety protec- 
tion. It can be obtained from pedigree and heterosis 
data, from morphological traits or using molecular 
markers which detect variation at the DNA sequence 
level (Smith and Smith 1992). In particular, DNA- 
based polymorphisms are a powerful tool in the assess- 
ment of the genetic similarity between breeding stocks 
(reviewed in Lee 1995). 

The discrimination power of restriction fragment 
length polymorphisms (RFLPs) has been extensively 
studied in maize, as has their use in establishing rela- 
tionships with yield and heterosis (Melchinger 1993). 
However, there are several drawbacks to RFLPs that 
have stimulated the development of alternative marker 
systems: large quantities of DNA are in fact required 
for RFLP analysis, which is costly, and the technique is 
difficult to automate. Moreover, it requires sizeable 
laboratories and specialised equipment. 

Various PCR-based marker techniques have recently 
been successfully introduced in the fingerprinting of 
plant genomes (Welsh and McClelland 1990; Kesseli 



et al. 1994) and in genetic diversity studies (Tinker et al. 
1993). Among them, random amplified polymorphic 
DNA (RAPD) analysis is quick (Welsh and McClelland 
1990; Williams et al. 1990) and well adapted for the 
efficient non-radioactive DNA fingerprinting of geno- 
types (dos Santos et al. 1994; Thormann et al. 1994). 
Problems with the reproducibility of amplification and 
with the scoring of error data have been reported for 
RAPDs (Demeke et al. 1997; Karp et al. 1997). 

Eukaryotic genomes are interspersed with tandem 
repeats of DNA, referred to as microsatellites or simple 
sequence repeats (SSRs). SSR polymorphisms have 
been extensively used as genetic markers in mammals 
(Tautz 1989); they occur frequently also in plant 
genomes, showing an extensive variation in different 
individuals and accessions (Akkaya et al. 1992; Senior 
and Heun 1993; Wu and Tanksley 1993). SSR loci are 
co-dominant markers more informative than RAPDs 
and RFLPs (Russell et al. 1997). Specific technical de- 
velopments are underway (Mitchell et al. 1997) that 
should result in the provision of SSRs that will be 
faster, more standardised and more effective than 
RFLP technology. 

Amplified fragment length polymorphism (AFLP™) 
is a multilocus marker technique developed by Vos 
etal. (1995). AFLP markers are genomic fragments 
detected after selective PCR amplification which pro- 
vide a number of appealing features in the fingerprint- 
ing of genomes of different complexity, including that 
of maize (Vos et al. 1995). The AFLP technique has 
been used to identify markers linked to disease resist- 
ance loci (Becker et al. 1995; Cervera et al. 1996), to 
fingerprint DNAs (Vos et al. 1995; Sharma et al. 1996), 
and to assess relationships between molecular poly- 
morphism and hybrid performance in maize (Ajmone- 
Marsan et al. 1998). 

A comparison of different marker techniques is 
timely, even though the utility of different molecular 
markers for soybean and barley germplasm has already 
been reported (Powell et al. 1996; Russell et al. 1997). 
The objectives of the present study were: (1) to compare 
the informativeness of different molecular markers and 
their applicability for genetic diversity analysis, geno- 
type identification and variety protection purposes, 
(2) to determine the genetic similarity obtained with 
RFLP- and PCR-based techniques in a set of maize 
inbred lines, and (3) to compare their effectiveness in 
estimating genetic similarity among maize inbreds. 



Materials and methods 

Plant materials and DNA extraction 

Thirty three inbred lines were chosen to explore the diversity of 
maize germplasm. All these inbreds have been extensively used in the 
production of hybrid seed and in maize breeding programs. Pedigree 
information was previously described in Livini et al. (1992). Based on 
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available information and on the heterotic behaviour in crosses, 13 
(A641, B14 A, B37, B73, B84, Cml09 Lol016, Lo916, Lo950, Lo951, 
Lo964, Lo999, and N28) can be associated with the Iowa Stiff Stalk 
Synthetic (BSSS) heterotic group, 13 (A619, C103, C123, H99, 
L088I, Lo924, Lo976, Lol077, Mol7, Oh43, Va22, Va59, and Va85) 
with the Lancaster Sure Crop (LSC), two to Wf9 (WF9 and Pa91), 
three to W153R (W153R, Lo932, and Lo944), and two to HY (H55 
and H96). Genomic DNA was isolated from a bulk of 20-30 shoots 
of 7-9-day old germinated seedlings and extracted using the CTAB 
method as previously described (Livini et al. 1992). 



Nucleic-acid manipulation and molecular-marker assays 

Conditions for restriction enzyme digestion, gel electrophoresis for 
RFLP, Southern transfer hybridisation, and autoradiography fol- 
lowed Livini et al. (1992). Forty seven genomic clones from the 
UMC and BNL collections and two restriction enzymes (EcoKl and 
i/mdIII) were used to characterise 53 RFLP loci in the 33 inbred 
lines. A total of 253 RFLP bands were binary coded as 1 or 0 for the 
presence or absence of such loci in each line, respectively. 

RAPD amplification was performed as described by Ajmone- 
Marsan et al. (1993) using a Perkin Elmer 9600 Thermal Cycler. 
Reaction products were analysed by electrophoresis on a 2% 
agarose gel and stained with ethidium bromide. A total of 85 primers 
(Operon Technologies, California, USA) were surveyed in two in- 
bred lines. Twenty five primers showing reproducible, and clearly 
scorable, polymorphic (present or absent) fragments, were used to 
fingerprint the 33 inbreds. 

The primers for the SSR markers were synthesised according to 
the sequences published by Senior and Heun (1993) and Taramino 
and Tingey (1996). SSR procedures were those described by 
Taramino and Tingey (1996). Processed fragments, along with load- 
ing dye and internal size standards, were run out on a 6% acrylam- 
ide gel (Pfeiffer et al. 1997) using an Automated Laser Fluorescent 
sequencing electrophoresis unit (Pharmacia). Fragments were label- 
led with fluorescein by direct incorporation of F-12-dUTP (2-uM 
final concentration) during the PCR reaction. Un-incorporated 
labelled nucleotide was removed by ethanol precipitation prior to 
loading samples on the gel. Data were processed using Fragment 
Manager Software v. 1.1 (Pharmacia). SSR bands were sized first 
and then binary coded by 1 or 0 for their presence or absence in each 
line. 

AFLP marker analysis was according to Vos et al. (1995). Briefly, 
total genomic DNA (400 ng) was restricted with 5 U of £coRI (rare 
cutter) and 5 U of Msel (frequent cutter) (Pharmacia), and double- 
stranded adapters ligated to the fragment ends. The structure of 
the adapter sequences, pre-amplification, amplification and poly- 
acrylamide-gel electrophoresis conditions were as in Ajmone-Mar- 
san et al. (1998). Polymorphic amplification products were visualised 
by autoradiography and scored manually. All AFLP polymor- 
phisms were scored as dominant markers. 

All names of the RFLP probes and the nucleotide sequences of the 
primers used for the amplification of AFLP, RAPD, and SSR 
markers are available on request. 



Data analysis 

The average number of alleles per locus, the allele frequency, the 
expected heterozygosity (He), and the effective number of alleles per 
locus were calculated as reported by Morgante et al. (1994). The 
total number of effective alleles (Ne) surveyed by RFLP, RAPD, 
SSR, and AFLP analyses was calculated by summing the number of 
effective alleles of all the analysed loci as Ne = £ne (i). To compare 
the efficiency among the four methods, where RFLPs and SSRs 
generally detect multiple alleles and one band per assay, whereas 
RAPDs and AFLPs detect two alleles and multiple bands per assay, 
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an assay efficiency index (Ai) was calculated. Ai combines the effec- 
tive number of alleles identified per locus and the number of the 
polymorphic bands detected in each assay as Ai = Ne/P, where Ne 
is the total number of effective alleles detected and P is the total 
number of assays performed for their detection. 

The genetic similarities (GSs) from RFLP, RAPD, SSR, and 
AFLP data were calculated among all possible pairs of lines using 
the Dice similarity index as in Nei and Li (1979). The co-ancestry 
coefficient, / between lines related by pedigree, was calculated as 
previously reported (Ajmone-Marsan et al. 1992). Cluster analyses 
were based on similarity matrices obtained with the unweighted pair 
group method using arithmetic averages (UPGMA) (Rohlf 1990) 
and relationships between inbred lines were visualised as dendro- 
grams. For each dendrogram the co-phenetic coefficients between 
the matrix of genetic similarities and the matrix of co-phenetic 
values were computed using appropriate routines of the NTSYS-pc 
package. The significance of the co-phenetic correlation observed 
was tested using the Mantel matrix correspondence test (Mantel 
1967). 

The bootstrap procedure was employed to determine the samp- 
ling variance of the genetic similarities calculated from the data sets 
obtained with the different marker systems. All data, irrespective of 
the dual or multiallelic nature of the marker system, were scored in 
the form of a binary matrix. For each pair of inbreds, the Dice 
similarity index (GS) was calculated from the 2000 random sub- 
samples at different sample sizes (10, 50, 100, 150, 200, and all bands 
when the total exceeded 200). Bootstrap standard deviation esti- 
mates were based on 2000 samples. The calculations were performed 
with the SAS macro "BOOT* (Jackknife and Bootstrap Analyses, 
SAS Institute Inc.). 



Results 

Levels of polymorphisms 

The 33 inbred lines were surveyed with the four differ- 
ent marker systems. All of the molecular markers were 
able to uniquely fingerprint each of the inbred lines. 
The levels of polymorphism detected with each marker 
system and the index comparing their informativeness 
are reported in Table 1. The total number of assays 
ranged from only six primer combinations for AFLPs 
to 53 probe/enzyme combinations for RFLPs. The to- 
tal number of polymorphic bands identified ranged 
from 90 for RAPDs to 253 for RFLPs. An average 
number of 4.8 alleles per locus, with an average effective 



number of 3.2 alleles per locus, ranging from 1.2 to 6.5, 
could be distinguished for each probe/enzyme combi- 
nation using RFLPs. This value increased to 6.8 with 
SSRs, with an average number of effective alleles of 4.4 
per locus, ranging from 1.1 to 6.6, while for RAPDs and 
AFLPs these values were lower (1.6 for both). This was 
reflected also in lower expected heterozygosity values. 
Overall the highest assay efficiency index was observed 
for AFLPs (61.9) and the lowest for RFLPs (3.2). 
RAPDs and SSRs (5.8 and 4.4, respectively) were com- 
parable to RFLPs. In particular, for AFLPs the high 
assay efficiency index is due to the simultaneous detec- 
tion of several polymorphic bands in a multiplex ampli- 
fication per single reaction. 



Genetic similarity 

A summary of the genetic similarity estimates, cal- 
culated for each marker system, between pairs of lines 
of the various heterotic groups is shown in Table 2. All 
marker systems indicated that lines of BSSS origin were 
more similar in comparison to inbred lines of other 
heterotic groups. The mean value of the GS estimate 
was, as expected, lower for BSSS x LSC crosses than 
within the BSSS and LSC groups themselves. This is 
consistent with the common practice in maize breeding 
of preferentially developing hybrids between heterotic 
groups because they are expected to perform better 
than those from crosses within heterotic groups. The 
similarity ranged from 0.92 within LSC types, using 
AFLPs, to 0.00 within BSSS x miscellaneous types, us- 
ing SSRs. The estimates of GS follow the same pattern 
across marker systems, i.e. higher estimates of sim- 
ilarity within the BSSS types and lower estimates 
within the LSC x miscellaneous types. Overall SSRs 
revealed the lowest similarity values and AFLPs the 
highest. 

The genetic similarity trees produced from each 
marker system are presented in Fig. 1. In these trees 
inbreds were ordered as expected, though with 
exceptions, into the major groups BSSS and LSC. 



Table 1 Level of polymorphism and comparison of informativeness obtained with RFLP, RAPD, SSR and AFLP markers in 33 maize 
inbred lines 



Parameters 


Marker system 










RFLP 


RAPD 


SSR 


AFLP 


Number of assay units 
Number of polymorphic bands 
Number of loci 

Average number of alleles per locus 
Expected heterozygosity 
Effective number of alleles per locus 
Assay efficiency index 


53 (probe/enzymes) 
253 
53 

4.8 

0.63 

3.2 

3.2 


25 (primers) 

90 

90* 

2.0 

0.36 

1.6 

5.8 


27 (primer pairs) 
183 
27 

6.8 

0.72 

4.4 

4.4 


6 (primer combination) 

232 

232 a 

2.0 

0.34 

1.6 
61.9 



"Theoretical maximum number of loci 



1251 

Table 2 Mean, minimum and maximum of the Dice genetic similarity coefficient (GS) calculated from different molecular marker systems for 
various groups of maize inbred lines 



Marker system GS value 







Within heterotic groups 




Between heterotic groups 




Between 
all lines 
ALL 
(n = 528) 






BSSS 
(n = 78) 


LSC 
(n = 78) 


MISC 
(»-21) 


BSSS x LSC 
(n = 169) 


BSSS x MISC 
(n = 91) 


LSC x MISC 
(ii-91) 


RFLP 


Min. 
Max. 
Mean 


0.28 
0.79 
0.48 


0.25 
0.86 
0.40 


0.90 
0.41 


0.18 
0.59 
0.36 


0.20 
0.49 
0.36 


fl 1 ft 

U. 18 

0.56 
0.35 


U.lo 

0.90 
0.37 


RAPD 


Min. 
Max. 
Mean 


0.34 
0.90 
0.64 


0.38 
0.80 
0.57 


0.31 
0.89 
0.53 


0.37 
0.76 
0.58 


0.36 
0.65 
0.53 


0.29 
0.68 
0.53 


0.29 
0.91 
0.56 


SSR 


Min. 
Max. 
Mean 


0.13 
0.84 
0.38 


0.08 
0.88 
0.31 


0.08 
0.82 
0.28 


0.07 
0.52 
0.27 


0.00 
0.52 
0.22 


0.04 
0.40 
0.20 


0.00 
0.88 
0.26 


AFLP 


Min. 
Max. 
Mean 


0.45 
0.89 
0.62 


0.43 
0.92 
0.59 


0.43 
0.89 
0.55 


0.36 
0.61 
0.48 


0.43 
0.64 
0.52 


0.39 
0.62 
0.51 


0.36 
0.92 
0.53 



Discrepancies in forming subgroups within the major 
groups were observed as well as in the clustering of 
inbred lines of miscellaneous origins. Considering the 
BSSS-related lines the topology of each tree is unique 
with some evident similarity: the clustering of B14, B37, 
and B73 is, for example, fully conserved. On the LSC 
side, clustering was consistently reported by all 
methods with the exception of Va22, a line derived 
from C103, indicating that all methods aggregated lines 
of different origin. The Oh43-related lines (Oh43 and 
A619) were positioned within the Lancaster group only 
by RAPDs and AFLPs, while SSRs and RFLPs clus- 
tered these with BSSS lines (although Oh43 is usually 
considered a Lancaster type). Similarly to the RAPD- 
based tree, the clustering based on AFLP data produc- 
ed a tree with a relatively narrow range of similarity 
values between the more-related and the more-distant 
pairs of inbreds. In spite of this, all the main clusters 
were confirmed by AFLP data. 

Four pairs of very similar inbreds (B14 A, Cml09; 
Lo932, Lo944; A619, Oh43; and H55, H96) were clus- 
tered together by all marker systems, while in one 
additional case similar lines (Lo916 and Lo999) were 
consistently grouped in AFLP, SSR and RAPD trees. 



Comparison between marker systems 

All the estimates of correlation coefficients (rs) among 
available co-ancestry coefficients (/ s) and genetic sim- 
ilarity (GS) data were highly significant (P < 0.01). 
RAPDs showed the lowest correlation (r = 0.40) with 
/ values, RFLPs and SSRs intermediate values 
(r = 0.57 and r = 0.53, respectively), while AFLPs 
showed the highest value (r = 0.62). The rs among 



similarity data were also significant. Correlation coeffi- 
cients of RAPD marker data (r = 0.51, r = 0.57, and 
r = 0.52 with RFLP, SSR and AFLP, respectively) with 
those obtained using other marker systems were lower 
than those among similarity estimates based on 
AFLPs, RFLPs, and SSRs (r = 0.70, r = 0.67 and 
r = 0.59, respectively between AFLP and RFLP, 
AFLP and SSR, and RFLP and SSR). The extent to 
which similarities were correlated varied considerably 
across the whole data set. When the set of pairwise data 
(528) was divided into two groups (according to the 
arithmetic mean of the observed GS range based on 
RFLP data: "more similar" lines with GS > 0.37 and 
"less similar" lines with GS < 0.37), it became evident 
that the genetic similarities estimated by different 
marker systems were mainly correlated only among 
similar lines, while the relationships among dissimilar 
lines were low and not significant. The GS values plot- 
ted against the estimate of co-ancestry value based on 
pedigree data followed the same pattern. 

The co-phenetic correlation coefficients provided for 
each marker system indicate the extent to which the 
clustering of genotypes depicted in the trees accurately 
represents the estimates of genetic similarity between 
inbreds obtained with that marker system. Overall the 
co-phenetic coefficients were medium to high, with the 
RFLP (0.84) and AFLP (0.83) data resulting in the 
highest correlation, SSR (0.80) showing an intermediate 
value, and the RAPD (0.72) assay producing the lowest 
correlation. 

All methods could clearly distinguish all 33 inbred 
lines, although the SSR data provided the highest level 
of discrimination between any pair of inbreds. In 
general, the grouping agreed with the pedigree informa- 
tion of the lines, although some discrepancies were 
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Fig. 2 Variation assessed by bootstrap sampling of genetic sim- 
ilarity between maize inbreds across different marker systems due to 
the different number of markers (bands) 



observed. In particular, genetic similarities based on 
AFLP data had the highest correlation with pedigree 
data, while those based on RAPDs had the lowest one. 



Bootstrap analysis 

To determine the sampling variance of genetic similar- 
ities calculated from different molecular marker data 
sets, bootstrap analysis with a declining number of 
bands was performed. The relationships between num- 
ber of bands and the sampling variance of the genetic 
similarity among all pairs of inbred lines for each 
method is presented in Fig. 2. The results indicated that 
above 150 bands there was a diminishing return in the 
precision gained by adding additional bands. As the 
number of bands moves below these thresholds the 
standard deviation begins to increase (and precision 
decreases) at a greater rate. 



Discussion 

In this paper we have shown that the number of alleles 
detectable in maize by SSRs is higher in comparison to 
other methods. This high level of polymorphism is to be 
expected because of the unique mechanism responsible 
for generating SSR allelic diversity by replication slipp- 
age (Tautz et al. 1986). It is also known that when SSRs 
have been compared to other marker systems they have 
revealed the highest level of polymorphism (Wu and 
Tanskley 1993; Morgante et al. 1994; Powell et al. 
1996). The present data indicate that on average SSRs 
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carry two-fold more information than AFLPs and 
RAPDs, and 40% more information than RFLPs, 
when the number of alleles per locus is the target. 

In agreement with previous observations (Becker 
etal. 1995), the lowest degree of polymorphism was 
associated with AFLPs. Conversely, the information 
measured as the assay efficiency index, which correlates 
with the number of effective alleles identified per assay, 
was more than ten-fold higher for AFLPs compared to 
the other methods. These findings are in good agree- 
ment with previous germplasm analysis carried out in 
several crop species (Lu et al. 1996; Powell et al. 1996). 
It can be concluded that SSRs are capable of revealing 
the highest level of information per single marker and 
that AFLPs detect the highest number of polymor- 
phisms in a single assay. This high assay efficiency 
index is a reflection of the efficiency of AFLPs to 
simultaneously analyse a large number of bands rather 
than the levels of polymorphism detected at each locus. 
The assay efficiency index for SSRs can, however, be 
considerably higher if multiplex PCR and gel-running 
procedures are adopted, where several microsatellites 
are simultaneously amplified and co-electrophoresed 
using multicolour fluorescent technologies (Lindqvist 
et al. 1996; Hey en et al. 1997). An additional advantage 
of SSRs over AFLPs is only relevant when mapping 
populations are derived from outcrossing heterozygous 
individuals, where the multi-allelism of SSR markers 
increases the number of informative genotypic classes 
in the progenies compared to the binary AFLP 
markers. 

The results shown by genetic similarity trees indicate 
that, except for RAPDs, they are highly similar. In 
addition, trees from these molecular methods agree 
with the information obtained from pedigree data. Sim- 
ilarly, Powell et al. (1996) found the lowest correlations 
among RAPDs and other marker systems. In this re- 
spect, it has been shown that RAPD analysis, based on 
the use of random primers, is likely to suffer from a lack 
of reproducibility due to mismatch annealing (Neale 
and Harry 1994). In the trees obtained from cluster 
analysis, all lines with defined affiliation to one of the 
heterotic groups were assigned to their specific main 
clusters, in agreement with the available data for maize 
(Ajmone-Marsan et al. 1992; Livini et al. 1992; Mumm 
and Dudley 1994; Smith et al. 1997). A second observa- 
tion is that, within the clusters, the grouping of more 
distantly related lines does not match precisely with the 
expectations based upon pedigree data. Differences 
among marker techniques in grouping genetically more 
distant lines have been previously reported (Powell 
et al. 1996). Other studies in Brassica and from pea 
accessions show that molecular marker-based similarit- 
ies and trees were significantly correlated across a wide 
range of germplasms (Thormann et al. 1994; Lu et al. 
1996). Many potential reasons for these discrepancies 
exist, including underlying assumptions in calculating 
pedigree data (Messmer et al. 1993), genome sampling 
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(Nei 1987), and the numbers of markers or probes 
employed (Tivang et al. 1994). 

The number of loci required for a reliable estimate of 
genetic similarity has been shown to vary from 15 
RFLP probes, giving 56 bands in Brassica sp. (dos 
Santos et al. 1994), to 100 RFLP clone-enzyme combi- 
nations (Messmer et al. 1993). Similarly, Tivang et al. 
(1994), investigating in maize the sampling variance of 
a RFLP data set in maize, found that the number of 
bands required for a CV of 10% was 388, 150, and 38 
for closely, intermediately, and distantly related in- 
breds, respectively. Our results using the bootstrap 
procedure suggest that 150 bands are sufficient for 
reliable estimates of genetic similarity. Accordingly, the 
average number of assays that could have been used in 
this study to attain such a precision in the estimate 
were 30-40 clone-enzyme combinations for RFLPs, 
40-50 primers for RAPDs, 20-30 primers for SSRs, 
and 4-5 enzyme combinations for AFLPs. Based on 
these estimates, the disagreement of the RAPD results 
in comparison to the other types of markers might be 
explained by the insufficient number of primers used. 

In conclusion, the results of this study indicate that, 
with the exception of RAPDs, the other DNA markers 
provide consistent information for germplasm identi- 
fication and pedigree validation. We have shown that 
SSR and AFLP profiling technologies can be good 
candidates to replace RFLP markers in genetic sim- 
ilarity estimates and variety description, and that they 
have comparable accuracy in grouping inbred lines 
selected by pedigree. They are generally much simpler 
to apply and more sensitive than the traditional 
morphological and biochemical methods or the 
RFLP-based fingerprinting techniques; yet they pro- 
vide results correlated with those from RFLP analyses. 
A major advantage of the SSR and AFLP methods is 
that they can be automated. While SSRs, thanks to 
their multi-allelism and co-dominance, appear to be 
suited for the analysis of outcrossing heterozygous indi- 
viduals, AFLPs, with their high multiplex ratio, offer 
a distinctive advantage when genome coverage is a ma- 
jor issue due to the presence of linkage disequilibrium, 
such as in inbred lines and breeding materials. 

Acknowledgements I. Pejic thanks the University of Udine for pro- 
viding his fellowship in 1994 and 1995. This work was supported by 
Ministero delle Risorse Agricole, Ambientali e Forestall, Roma, 
Italy, special grant: Piano Nazionale Biotecnologie Vegetali. 



References 

Ajmone-Marsan P, Livini C, Messmer MM, Melchinger AE, Motto 
M (1992) Cluster analysis of RFLP data from related maize 
inbred lines of the BSSS and LSC heterotic groups and compari- 
son with pedigree data. Euphytica 60: 139-148 

Ajmone-Marsan P, Egidy G, Monfredini G, Di Silvestro S, Motto 
M (1993) RAPD markers in maize genetic analysis. Maydica 
38:259-264 



Ajmone-Marsan P, Castiglioni P, Fusari F, Kuiper M, Motto 
M (1998) Genetic diversity and its relationship to hybrid perfor- 
mance in maize as revealed by RFLP and AFLP markers. Theor 
Appl Genet 98:219-227 

Akkaya MS, Bhagwat A A, Cregan PB (1992) Length polymorphism 
of simple sequence repeat DNA in soyabean. Genetics 
132:1131-1139 

Becker J, Vos P, Kuiper M, Salamini F, Heun M (1995) Combined 
mapping of AFLP and RFLP markers in barley. Mol Gen Genet 
249:65-73 

Cervera MT, Gusmao J, Steenackers M, Peleman J, Storme V, 
Venden Broeck A, Van Montagu M, Boerjan W (1996) Identi- 
fication of AFLP molecular markers for resistance against 
Melampsora larici-populina in Populus. Theor Appl Genet 93 : 
733-737 

Demeke T, Sasikumar B, Hucl P, Chibbar RN (1997) Random 
amplified polymorphic DNA (RAPD) in cereal improvement. 
Maydica 42:133-142 

Hallauer AR, Russell WA, Lamkey KR (1988) Corn breeding. In: 
Sprague GF, Dudley JW (eds) Corn and corn improvement, 
3rd edn. Agron Monogr 18. ASA, CSS A, and SSSA, Madison, 
Wisconsin, USA 

Heyen DW, Beever JE, Da Y, Evert RE, Green C, Bates SR, Ziegle 
JS, Lewin HA (1997) Exclusion probabilities of 22 bovine micro- 
satellite markers in fluorescent multiplexes for semiautomated 
parentage testing. Anim Genet 28 : 21-27 

Karp A, Edwards K, Bruford M, Vosman B, Morgante M, Seberg O, 
Kremer A, Boursot P, Arctander P, Tautz D, Hewitt G 
(1997) Newer molecular technologies for biodiversity evalu- 
ation: opportunities and challenges. Nature Biotechnol 
15:625-628 

Kesseli R, Ochoa O, Michelmore R (1994) Variation of RFLP loci in 
Lactuca spp. and origin of cultivated lettuce (L. sativa). Genome 
34:430-436 

Lee M (1995) DNA markers and plant breeding programs. Adv 
Agron 55:265-344 

Lindqvist AK, Magnusson PK, Balciuniene J, Wadelius C, Lind- 
holm E, Alarcon-Riquelme ME, Gyllensten UB (1996) Chromo- 
some-specific panels of tri- and tetra-nucleotide microsatellite 
markers for multiplex fluorescent detection and automated 
genotyping: evaluation of their utility in pathology and forensics. 
Genome Res 6 : 1 170-1 176 

Livini C, Ajmone-Marsan P, Melchinger AE, Messmer MM, Motto 
M (1992) Genetic diversity of maize inbred lines within and 
among heterotic groups revealed by RFLP. Theor Appl Genet 
84:17-25 

Lu J, Knox MR, Ambrose MJ, Brown JKM, Ellis THN (1996) 
Comparative analysis of genetic diversity in pea assessed by 
RFLP- and PCR-based methods. Theor Appl Genet 93: 
1103-1111 

Mantel N (1967) The detection of disease clustering and a generaliz- 
ed regression approch. Cancer Res 27 : 209-220 

Melchinger AE (1993) Use of RFLP markers for analysis of genetic 
relationship among breeding materials and prediction of hybrid 
performance. In: Buxton DR et al. (eds) International Crop 
Science I. CSS A, Madison, Wisconsin, USA, pp 621-628 

Messmer MM, Melchinger AE, Hermann RG, Boppenmeier J 
(1993) Relationship among early European maize inbreds. 
II. Comparison of pedigree and RFLP data. Crop Sci 
33:944-950 

Mitchell SE, Kresovich S, Jester CA, Hernandez CJ, Szewc-McFad- 
den AK (1997) Application of multiplex PCR and fluorescence- 
based semi-automated allele-sizing technology for genotyping 
plant genetic resources. Crop Sci 37:617-624 

Morgante M, Rafalsky A, Biddle P, Tingey S, Olivieri AM (1994) 
Genetic mapping and variability of seven soyabean simple se- 
quence repeat loci. Genome 37 : 763-769 

Mumm RH, Dudley JW (1994) A classification of 148 U.S. maize 
inbreds. I. Cluster analysis based on RFLPs. Crop Sci 34: 
842-851 



Neale DB, Harry DE (1994) Genetic mapping in forest trees: RFLPs, 

RAPDs and beyond. AgBiotech News Inf 6: 107N-114 N 
Nei M (1987) Molecular evolutionary genetics. Colombia University 

Press, New York 
Nei M, Li WH (1979) Mathematical model for studying genetic 

variation in terms of restriction endonucleases. Proc Natl Acad 

Sci USA 76:5269-5273 
Pfeiffer A., Olivieri AM, Morgante M (1997) Identification and 

characterization of simple sequence repeats in Norway spruce 

(Picea abies K.). Genome 40:411-419 
Powell W, Morgante M, Andre C, Hanafey M, Vogel J, Tingey S, 

Rafalsky A (1996) The comparison of RFLP, RAPD, AFLP and 

SSR (microsatellite) markers for germplasm analysis. Mol Breed 

2:225-238 

Rohlf FJ (1990) NTSYS-pc numerical taxonomy and multivari- 
ate analysis system, version 1.50. Exeter Publications, 
New York 

Russell JR, Fuller JD, Macaulay M, Hatz BG, Jahoor A, Powell W, 
Waugh R (1997) Direct comparison of levels of genetic variation 
among barley accessions detected by RFLPs, AFLPs, SSRs and 
RAPDs. Theor Appl Genet 95: 714-722 

Santos JB dos, Nienhuis J, Skroch P, Tivang J, Slocum MK (1994) 
Comparison of RAPD and RFLP genetic markers in determin- 
ing genetic similarity among Brassica oleracea L. genotypes. 
Theor Appl Genet 87 : 909-9 1 5 

Senior ML, Heun M (1993) Mapping maize microsatellites and 
polymerase-chain-reaction confirmation of the targeted repeats 
using a CT primer. Genome 36:884-889 

Sharma SK, Knox MR, Ellis THN (1996) AFLP analysis of the 
diversity and phylogeny of Lens and its comparison with RAPD 
analysis. Theor Appl Genet 93 : 751-758 

Smith JSC, Smith OS (1992) Fingerprinting crop varieties. Adv 
Agron 47:85-129 



1255 

Smith JSC, Chin ECL, Shu H, Smith OS, Wall SJ, Senior ML, 
Mitchell SE, Kresovich S, Ziegle J (1997) An evaluation of the 
utility of SSR loci as molecular markers in maize (Zea mays L.): 
comparison with data from RFLPs and pedigree. Theor Appl 
Genet 95:163-173 

Taramino G, Tingey S (1996) Simple sequence repeats for ger- 
mplasm analysis and mapping in maize. Genome 39 : 277-287 

Tautz D (1989) Hypervariability of simple sequences as a general 
source for polymorphic DNA markers. Nucleic Acids Res 
17:6463-6471 

Tautz D, Trick M, Dover GA (1986) Cryptic semplicity in DNA is 
a major source of variation. Nature 322 : 652-656 

Thormann CE, Ferreira ME, Camargo LEA Tivang JG, Osborn TC 
(1994) Comparison of RFLP and RAPD markers to estimate 
genetic relationship within and among cruciferous species. Theor 
Appl Genet 88:973-980 

Tinker NA, Fortin MG, Mather DE (1993) Random amplified 
polymorphic DNA and pedigree relationship in spring barley. 
Theor Appl Genet 85:976-984 

Tivang JG, Nienhuis J, Smith OS (1994) Estimation of sampling 
variance of molecular-marker data using the bootstrap proced- 
ure. Theor Appl Genet 89 : 259-264 

Vos P, Hogers R, Bleeker M, Reijans M, Van De Lee T, Homes M, 
Frijters A, Pot J, Peleman J, Kuiper M, Zabeau M (1995) AFLP: 
a new technique for DNA fingerprinting. Nucleic Acids Res 
23:4407-4414 

Welsh J, McClelland M (1990) Fingerprinting genomes using PCR 
with arbitrary primers. Nucleic Acids Res 18:7213-7218 

Williams JGK, Kubelik AR, Livak KJ, Rafalsky JA, Tingey SV 
(1990) DNA polymorphisms amplified by arbitrary primers are 
useful as genetic markers. Nucleic Acids Res 18 : 6531-6535 

Wu K-s, Tanksley SD (1993) Abundance, polymorphism and genetic 
mapping of microsatellites in rice. Mol Gen Genet 241 : 225-235 



Theor Appl Genet (1997) 95 : 714-722 © Springer- Verlag 1997 



J. R. Russell • J. D. Fuller • M. Macaulay 

B. G. Hatz • A. Jahoor • W. Powell • R. Waugh 



Direct comparison of levels of genetic variation among barley accessions 
detected by RFLPs, AFLPs, SSRs and RAPDs 



Abstract RFLPs, AFLPs, RAPDs and SSRs were used 
to determine the genetic relationships among 18 culti- 
vated barley accessions and the results compared to 
pedigree relationships where these were available. All 
of the approaches were able to uniquely fingerprint 
each of the accessions. The four assays differed in the 
amount of polymorphism detected. For example, all 13 
SSR primers were polymorphic, with an average of 5.7 
alleles per primer set, while nearly 54% of the frag- 
ments generated using AFLPs were monomorphic. The 
highest diversity index was observed for AFLPs (0.937) 
and the lowest for RFLP (0.322). Principal co-ordinate 
analysis (PCoA) clearly separated the spring types from 
the winter types using RFLP and AFLP data with the 
two-row winter types forming an intermediate group. 
Only a small group of spring types clustered together 
using SSR data with the two-row and six-row winter 
varieties more widely dispersed. Direct comparisons 
between genetic similarity (GS) estimates revealed by 
each of the assays were measured by a number of 
approaches. Spearman rank correlation ranked over 
70% of the pairwise comparisons between AFLPs 
and RFLPs in the same order. SSRs had the lowest 
values when compared to the other three assays. These 
results are discussed in terms of the choice of appropri- 
ate technology for different aspects of germplasm 
evaluation. 
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Introduction 

As one of the first crop plants to be domesticated, 
barley (Hordeum vulgare L.) remains one of the most 
important crops today. Ranking fourth in world acre- 
age, barley is used for human consumption, as a fodder 
crop and as a raw material for brewing beer and whisky 
(Brown 1992). It belongs to the genus Hordeum, which 
comprises over 32 species, including diploid and poly- 
ploid, perennial and annual types, which are spread 
throughout the world. The genus can be divided into 
three groups of varying importance to cultivated barley 
improvement; the primary genepool (H. vulgare spp. 
vulgare and H. vulgare spp. spontaneum), the secondary 
genepool (if. bulbosum) and the tertiary genepool (all 
other Hordeum species). Presently more than 250,000 
Hordeum accessions are held in genebanks through- 
out the world, and the number is increasing (IBPGR 
1992). With the growth of the germplasm collection 
a need for procedures which will allow their more 
effective use is required. The 'Core Collection' concept, 
is one such method, which should provide users with 
a limited set of genetically distinct and representative 
accessions (Brown 1989). Recently this concept has 
been applied to barley, and the Barley Core Collection 
(BCC), consisting of a limited sample of accessions 
considered to represent the spectrum of genetic diver- 
sity available in the genus, was established (Hintum 
1992). 

In such collections, morphological data are the prin- 
ciple descriptors which have been used to detail the 
accessions held. With the development of molecular 
markers and their many perceived advantages, it is 
crucial that these techniques are applied to assess 
genetic diversity in germplasm collections in order 
to supplement and refine the morphological-based 
classification. However, in recent years, the number 
of molecular assays available for application in this 
area has increased dramatically, with each method dif- 
fering in principle, in application, in the type and 
amount of polymorphism detected and in cost and time 
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requirements. The approaches include restriction frag- 
ment length polymorphism (RFLPs; Botstein et al. 
1980), random amplified polymorphic DNA (RAPDs; 
Williams et al. 1990), simple sequence repeat polymor- 
phisms or microsatellites (SSRs; Tautz 1989) and Am- 
plified Fragment Length Polymorphism (AFLPs; 
Zabeau and Vos 1993). 

Faced with this wealth of marker technology, it is 
appropriate to determine if the same patterns of varia- 
bility are revealed by each and whether the observed 
molecular diversity reflects either co-ancestry or mor- 
phological classification. To address this we have 
evaluated and compared similarity measures obtained 
from the four above systems on a set of accessions 
which are representative of cultivated European barley 
germplasm. This has allowed us to compare the results 
obtained from molecular analysis with each other and 
with pedigree information. The results are discussed in 
relation to the overall genetic diversity observed and 
the features of the individual assays. 



Marker analysis 
RAPD 

RAPD amplifications were performed as described by Barua et al. 
(1993). Fragments were separated on 1.5% agarose gels, stained with 
ethidium bromide, visualised with ultraviolet light and photo- 
graphed. The presence or absence of polymorphic bands were 
scored. Twenty primers, which were polymorphic between the par- 
ents of a spring x spring cross ('Blenheim* and E224/3), were used in 
this study. 



RFLP 

RFLP profiles were detected according to the protocol described by 
Graner et al. (1991). DNA was digested with three restriction en- 
zymes (BamHl, EcoKL and JfmdIII), and restriction fragments were 
detected using 48 single-copy DNA clones selected from previous 
mapping experiments to give good genome coverage and levels of 
polymorphism (Graner et al. 1991). RFLP patterns were scored as 
presence or absence of bands. 



Materials and methods 

Plant material and DNA isolation 

Eighteen accessions (Table 1), representing the majority of an- 
cestors European cultivated barley, were selected for this study. 
Total genomic DNA was isolated from fresh leaf material by a 
modification of the method described by Saghai-Maroof et al. 
(1984). 



SSR 

Two sources of simple sequence repeats were used in this study: 
database-derived repeats and repeats derived from an enriched 
genomic library. The 6 database-derived SSRs are described in 
a recent publication by Becker and Heun (1995). The 7 library- 
derived SSRs are described by Macaulay et al. (in preparation). SSR 
assays were performed as described by Morgante et al. (1994). Allele 
lengths were determined by comparing the most intense band with 
an M13 DNA sequence marker. 



Table 1 Pedigree information 
and country of origin of 18 barley 
accessions used in molecular 
analysis 



Cultivar 



Pedigree 



Origin* 



Spring type, two-rowed: 
Aramir 
Beka 

Golden Promise 
Grit 

Hora 
Krona 
Triumph 

Union 
Volga 

Winter type, two-rowed: 
Igri 

Marinka 
Romanze 
Sonja 

Winter type, six-rowed: 
Borwinia 
Express 
Franka 
Gaulois 
Rondo 



VollaxEmir NL 

Bethge XIII x Kniefel F 

X-ray mutant from Maythorpe (Irish Goldthorpe x Maja) GB 

Langenstein-Nungesser (5547/67 x 46459/68) 480/68 or Hadml. D 
554-Emir-11191-Union-46495-Diamant 14008 

Sultan x (Weihenstephaner 1206 Nacktgerste x Volla) D 

Complex cross including Triumph D 

(Hadm.24566 x Diamant x 1402964/6) x ((Alsa x Abyssinian) D 
x St. x Union) 

(Weihenstephaner Mehtauresistente II x Donaria) x Firlbecks III D 

Complex cross with eight varieties F 

(Malta x Carlsberg 1427) x Ingrid D 

(Alpha x SVP 674) x Malta NL 

Weihenstephan 4622/73 x (Malta x Sonja) D 

Tria x Malta D 

Vogelsanger Gold x St. 7246 D 

Robur x Athene D 

(Vogelsanger Gold x Senta) x (Dura x Dea) x Vogelsanger Gold D 

Gerbel x Athene F 

Tanaroo x Sisfor L. 90 I 



a NL, The Netherlands; F, France; GB, Great Britain; D, Germany; I, Italy 
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AFLP 

AFLP analysis was essentially as described by Vos et al. (1995). 
Briefly, 500 ng of genomic DNA was digested with EcoKl and Msel 
and double-stranded adaptors ligated to the fragment ends. This was 
followed by a pre-amplification step using non-selective primers. 
Selective amplifications were performed on the pre-amplified frag- 
ment mixture using a total of six primer combinations. Only the 
EcoRI primer was radiolabeled with y-[ 33 P] ATP (ICN), and all 
primers had three selective nucleotides. Amplification products were 
separated by denaturing 6% polyacrylamide gel electrophoresis 
(PAGE), visualised by autoradiography and manually scored for the 
presence or absence of bands. 

All of the primer names and sequences used are available on 
request from the authors. 



Data analysis 

Diversity values were calculated for each locus as (1 - 2 Pi 2 , where 
Pi is the phenotypic frequency for each assay unit (RFLPs- 
probe/enzyme combinations; RAPDs-primers; SSRs-primer pairs; 
AFLPs-primer combinations). Genetic similarities (GS) were cal- 
culated using the GENSTAT Version 531 software package accord- 
ing to Nei and Li*s (1979) estimate of similarity. Similarities were 
expressed using the group average agglomerative clustering function 
of GENSTAT to generate principal co-ordinate plots (Kempton and 
McNicol 1990). Correlations between assays were calculated using 
[Procrustes rotational analysis (PR) on the principal co-ordinate 
data] Spearman rank correlation (SRC) and linear regression of the 
GS values. 



Results 

Fingerprinting 

All of the molecular approaches used in this study were 
able to uniquely fingerprint each of the 18 cultivated 
barley accessions. The total number of assay units 
varied for each marker system from only 6 primer 
combinations for AFLPs to 144 probe/enzyme combi- 
nations for RFLPs (Table 2). Similarly, the number of 
bands scored ranged from 70 for SSRs to 299 for 
RFLPs. The percentage of polymorphic bands for each 
assay did not correlate to the total number of bands. 
For example, only 70 bands were scored for SSRs, 
which was the lowest number, but all 70 were polymor- 
phic. In contrast, 297 AFLP bands were scored, and 
only 46.8% of those were polymorphic. RFLPs and 



RAPDs were intermediate with 83.2% and 66.3%, re- 
spectively, of all bands scored being polymorphic. 
There was wide variation in the average number of 
genotypes revealed by each marker system (Fig. 1). 
With RFLPs, for each probe/enzyme combination, 
an average of 2.37 genotypic classes could be dis- 
tinguished. With AFLPs this figure increased to 17.2 
as nearly all primer combinations were able to dis- 
criminate between the 18 accessions used. This is 
further reflected in the diversity index measures. 
Overall the highest diversity index was observed for 
AFLPs (0.937), and the lowest for RFLPs (0.322). 
RAPDs and SSRs were intermediate (0.521 and 0.566, 
respectively). 



Genetic similarity 

The cultivated barley genepool can be divided into 
spring and winter types. The winter barleys are mainly 
used for fodder and can be further divided into two- 
and six-rowed types. The spring barleys are mainly 
used for malting. The maximum, minimum and mean 
similarity estimates between the spring barleys and 
two-row and six-row winter barleys for each assay 
system are shown Table 3. The similarities ranged from 
0.97 within spring types using AFLPs to 0.45 within 
six-row winter types using SSRs. Between assay sys- 
tems the estimates of similarity followed the same pat- 
tern, i.e. higher estimates of similarity within the spring 
types (means: RFLPs = 0.843, AFLPs = 0.924, 
SSRs = 0.829) and lower estimates within the six-row 
winter types (means: RFLPs = 0.70, AFLPs = 0.877, 
SSRs = 0.657). Estimates with two-row winter types 
were intermediate. The situation with RAPDs was 
different, with spring and six-row winter types exhi- 
biting equivalent mean similarities (0.879 and 0.897, 
respectively). Overall, SSRs revealed the lowest sim- 
ilarity values (0.93-0.45) and AFLPs the highest 
(0.97-0.81). 

Some accessions can be traced to common ancestors. 
For example, Grit and Triumph have Union in their 
pedigrees and Krona has Triumph. Our expectation 
would therefore be that these 4 accessions should be 
closely related. Table 4 shows the genetic similarity 



Table 2 Analysis of the RFLP-, 
RAPD-, SSRs- and 
AFLP-generated banding 
patterns 



Marker 


Number of 


Total 


Number of 


Number of 


Number of 


Diversity 




assay units 


no of 


polymorphic 


bands 


phenotypes 


index 




bands 


bands (%) 


per assay 


per assay 












unit 


unit 




RFLPs 


114 (42 probes, 


299 


249(83.2%) 


2.62 


2.37 


0.322 




3 enzymes) 










0.521 


RAPDs 


22 (primers) 


107 


71 (66.3%) 


4.86 


3.41 


SSRs 


13 (primer pairs) 


70 


70 (100%) 


5.38 


5.38 


0.566 


AFLPs 


6 (primer 


297 


139(46.8%) 


49.5 


17.2 


0.937 




combinations) 
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Fig. 1A-C An example of the different information content ob- 
served with RAPDS (A), SSRs (B) and AFLPs (Q 

values for the comparisons of these 4 accessions with 
each of the molecular assays. With RFLPs, AFLPs and 
SSRs the genetic similarity values were higher than the 
mean values for all the spring types, and the highest 
similarity was between Triumph and Grit (RFLPs 0.93, 
AFLPs 0.97, SSRs 0.97). From the pedigree informa- 
tion in Table 1, Triumph and Grit share a number of 



parental lines including Union, Diamont and Hadm. 
With RAPDs, the genetic similarity values were less 
than the average, although the Triumph and Grit com- 
parison was again the highest. 

The genetic similarity values for the two-row winter 
varieties were intermediate between the spring and 
six-row winter types for RFLPs, AFLPs and SSRs. 
Sonja and Romanze were more similar than the other 
two-row winter types with values of 0.84 (SSRs), 0.93 
(AFLPs), 0.89 (RFLPs) and 0.94 (RAPDs). This was not 



Table 3 Maximum, minimum and mean genetic similarity estimates calculated from RFLP, RAPD, SSRs and AFLP data for winter and 
spring types 

RFLPs RAPDs AFLPs SSRs Parentage 

Max Min Mean Max Min Mean Max Min Mean Max Min Mean Max Min Mean 



Spring 


93.0 


76.0 


84.3 


95.0 


84.0 


87.9 


97.0 


Two-row 


89.0 


81.0 


83.8 


95.0 


86.0 


91.7 


93.0 


winter 
















Six-row 


85.0 


60,0 


70.0 


95.0 


88.0 


89.7 


91.0 



winter 



88.0 


92.4 


93.0 


66.0 


82.9 


0.330 


0.020 


0.133 


98.0 


91.0 


84.0 


56.0 


71.3 


0.290 


0.100 


0.212 


81.0 


87.7 


90.0 


45.0 


65.7 


0.350 


0.000 


0.111 
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Table 4 Genetic similarity values 
for the comparisons of 4 spring 
accessions with each of the 
molecular assays 





Grit 


Krona 


Triumph 


Union 


Maximum 


Minimum 


Mean 


T> A DHi- 
















Ont 


1 on rv 














K.rona 


OO. u 


100.0 












Triumph 


90.0 


85X) 


100.0 










Union 


88.0 


88.0 


84.0 


100.0 


95.0 


84.0 


87.9 


J> PI Do- 
Kr LrS. 
















Grit 


luu.u 














Krona 


r^ n 


1UU.U 












Triumph 


93.0 


88.0 


100.0 










Union 


87.0 


87.0 


86.0 


100.0 


93.0 


76.0 


84.3 


ArLrs: 
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lUU.U 














Krona . 


01 O 


i ah n 












Triumph 


97.0 


92.0 


100.0 










Union 


93.0 


94.0 


94.0 


100.0 


97.0 


88.0 


92.4 


SSRs: 
















Grit 


100.0 














Krona 


84.0 


100.0 












Triumph 


97.0 


87.0 


100.0 










Union 


69.0 


70.0 


94.0 


100.0 


97.0 


66.0 


82.9 



unexpected as the co-efficient of parentage values were 
also the highest (0.290 for Sonjax Romanze compared to 
the mean for two-row winter type of 0.212). Both Sonja 
and Romanze are related through Malta, and Romanze 
has Sonja in its pedigree. The lowest genetic similarity 
values were observed for comparisons with Rondo. 



Genetic relatedness 

Associations among the 18 accessions were revealed by 
principal co-ordinate analysis (PCoA) (Fig. 2). The 
PCoA for the combined data (775 bands) clearly separ- 
ated the winter from the spring accessions. Among the 
winter types, the two-rowed and six-rowed varieties 
formed two distinct groups, with the two-rowed types 
forming an intermediate group between the spring and 
six-rowed winter types. In the PCoAs generated by 
RFLP (299 bands) and AFLP (297 bands) data, a sim- 
ilar arrangement was observed. From the RAPD data, 
three distinct groups were again observed, although the 
spring types were more dispersed. Only a small group 
of spring types clustered together using SSR data, and 
two-row and six-row winter types were again more 
dispersed. On all of the PCoAs, Rondo appears in 
a remote position. In addition, 'Volga', a spring variety, 
was positioned between the rest of the spring and the 
two-rowed winter types. 



Comparison between assays 

To compare the results obtained with the four tech- 
niques, we tested correlations using Procrustes rotation 



(PR), linear regression of the pairwise GS values (LR) 
and Spearman rank correlation (SRC). The results for 
SRC (which compares how each system ranks pairwise 
similarities) are shown in Table 5. Comparisons using 
PR and LR showed the same general trends although 
the overall correlations were lower. Over 70% of the 
pairs of genotypes were ranked in the same order with 
RFLPs and AFLPs. This correlation is reduced to 
10.9% when comparing RAPDs with AFLPs. SSRs 
were intermediate with over 50% of the genotypes 
ranking in the same order as that obtained with AFLPs 
and RFLPs. 



Discussion 

Given the proliferation of genetic markers, compari- 
sons between techniques are inevitable. However, there 
is a need for such comparisons in order to decide on 
which technique is best suited to the issues being exam- 
ined. In this study, three of the newer polymerase chain 
reaction (PCR)-based systems (RAPDs, SSRs and 
AFLPs) developed during the last 5 years have been 
compared with the well established RFLP system that 
was developed over 15 years ago. Each technique not 
only differs in principal, but also in the type and 
amount of polymorphism detected. The levels of poly- 
morphism between the four techniques varied widely, 
ranging from a maximum of 100% (SSRs) to only 
48.6% (AFLPs). Similar results were observed when 
Rus-Kortekaas et al. (1994) directly compared SSRs 
with RAPDs in tomato where the level of polymor- 
phism was 40% with RAPDs compared to 100% with 
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Fig. 2 Associations among the springs and winters cultivars re- 
vealed by principal co-ordinate analysis for each molecular assay 

Table 5 Correlations obtained using RFLPs, RAPDs, AFLPs and 
SSRs based on Spearman's rank correlation and Procrustes rotation 
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SSRs. Indeed, whenever SSRs have been compared to 
other systems, they have always revealed the highest 
levels of polymorphism (Rus-Kortekaas et al. 1994; 
Salimath et al. 1995; Saghai Maroof et al. 1994; Powell 
et al. 1996; Maughan et al. 1995; Morgante et al. 1994; 
Wu and Tanksley 1993). The level of polymorphism 
detected using RFLPs in this study was higher (83.2%) 
than that observed in previous studies on barley using 
a similar selection of genotypes (46%) (Melchinger 
et al. 1994). This is probably due to pre-selection of 



polymorphic RFLP probes. The lowest level of poly- 
morphism was associated with AFLPs. Becker et al. 
(1995) also observed that levels of polymorphism re- 
vealed by AFLPs were lower than by RFLPs. How- 
ever, although AFLPs do not offer the highest level of 
polymorphism, they are the most efficient because they 
have the capacity to reveal many polymorphic bands in 
a single lane. The average number of bands per lane or 
per PCR for AFLPs was 49.5, compared to 1.0 band 
per lane or PCR for SSRs. Thus, when the overall 
diversity indices of the four techniques were compared, 
AFLP was the highest (0.937). Powell et al. (1996) 
introduced the concept of Marker Index as an overall 
measure of marker efficiency, and they demonstrated 
that, in Glycine, AFLPs had the highest Marker Index 
compared to other available marker systems. The high 
Marker Index or diversity index is a reflection of the 
efficiency of AFLPs to simultaneously analyse a large 
number of bands rather than the levels of polymor- 
phism detected. 

Barley germplasm can be divided into two 
genepools, winter and spring, based on morphology 
distinctions. Melchinger et al. (1994) using RFLPs ob- 
served a clear separation between the spring and winter 
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types. In this study, similar results were observed using 
RFLPs, AFLPs, RAPDs and SSRs. Furthermore, Mel- 
chinger et al. (1994) noted that sub-groups were also 
apparent for accessions with similar pedigrees, such as 
the compact grouping of two-row winter types inter- 
mediate between the spring and six-row winter types. 
With the exception of the SSR data, the two-row winter 
types form a sub-group between the six-row winter and 
the spring types for RFLP, AFLPs and RAPDs. With 
the SSR data there is a clear separation between the 
spring and winter types, but not within the winter 
types. This is not unexpected considering the low level 
of band sharing between accessions; even within groups 
the estimates of genetic similarity were much lower 
than any of the other assays. 

Several previous studies have compared the use of 
RFLPs and RAPDs to examine genetic relatedness 
(Hallden et al. 1994; Thormann et al. 1994; Liu and 
Furnier 1993; dos Santos et al. 1994), and most of these 
show that RAPDs and RFLPs detect very similar rela- 
tionships among the same group of accessions. Re- 
cently, other reports have compared RAPDs or RFLPs 
and SSRs on the same set of genotypes (Rus-Kortekaas 
et al. 1994; Wu and Tanksley 1993; Salimath et al. 1995; 
Maughan et al. 1995). Rus-Kortekaas et al. (1994) ob- 
served a lower percentage of band sharing in tomato 
accessions with SSRs compared to RAPDs and sugges- 
ted that higher band sharing would make RAPDs more 
suitable for genetic relatedness studies. The results in 
this study would support the finding that SSRs may not 
be particularly well suited for pedigree relationship 
studies, although only a small number of SSRs were 
used. 

Knowledge of genetic variation and the genetic 
relationship between genotypes is an important consid- 
eration for efficient rationalisation and utilisation of 
germplasm resources. Furthermore, it is important for 
the optimal design of plant breeding programmes, in- 
fluencing the choice of genotypes to cross for the devel- 
opment of new populations. In barley, breeders have 
made crosses between highly selected genotypes with 
the result that the number of genotypes within the 
breeding genepool is very small. According to Graner 
et al. (1994) better knowledge and measures of genetic 
similarity of accessions could help to maintain genetic 
diversity. In the past, indirect estimates of similarity 
based on pedigree information have been widely used 
in many species including barley. Such estimates may 
not always reflect the true relationships between acces- 
sions (Graner et al. 1994). In this study we have used 
molecular markers to determine direct measures of 
genetic similarity between individuals. The estimates 
varied from 0.97 (AFLPs) to 0.45 (SSRs). Melchinger et 
al. (1994) reported GS values of 0.79 for unrelated 
barley pairs, based on RFLPs. The RFLP results re- 
ported in this paper were similar to these. Also, Tinker 
et al. (1993) observed GS values in a set of 27 North 
American barley cultivars using RAPDs which were 



similar to those found here (0.84-0.95). The values of 
GS based on SSRs in this present study are much lower 
than those based on RFLPs, AFLPs and RAPDs. Rus- 
Kortekaas et al. (1994) reported that the percentage of 
band sharing between tomato cultivars using SSRs was 
only 50.8% compared with 82.7% for RAPDs. 
Plaschke et al. (1995) observed even lower (0.31) esti- 
mates of genetic similarity when employing SSRs to 
examine wheat accessions and suggested that these low 
values are a reflection of the high information content 
provided by SSRs. 

Although we have shown that molecular approaches 
can be used to group barley cultivars into mor- 
phologically distinct groups, and also further into 
sub-groups which have a similar genetic background, 
we have not addressed the issue of concordance of 
molecular-based estimates of GS and co-ancestry. 
Graner et al. (1994) compared RFLP-based estimates 
of GS with co-ancestry for a set of 48 cultivars. A very 
weak correlation was reported; r s = 0.21 for winter 
and r s = 0.42 for spring types. Similarly, using protein- 
based gliadin markers Cox et al. (1985) observed a 
correlation of r s = 0.27. Both Graner et al. (1994) 
and Cox et al. (1985) agree that perhaps the reason 
for these poor correlations may be the high back- 
ground similarity found for unrelated accessions using 
molecular markers. When related cultivars were used 
to investigate correlations between RAPD-based 
estimates of GS and co-ancestry a moderate correlation 
of r s = 0.61 was observed between both measures 
(Tinker et al. 1993). Plaschke et al. (1995) observed 
similar results in wheat using SSR-based GS estimates 
and pedigree measures (r s = 0.55). Although we have 
only a limited set of co-ancestry measures for the acces- 
sions studied here, several conclusions can be drawn 
from the correlations between molecular estimates of 
GS and the co-efficient of parentage. For example, 
the co-efficient of parentage for Rondo was 0 for all of 
the pairwise comparisons, and with all molecular 
measures Rondo had the lowest GS value. The low-to- 
moderate correlations between molecular measures of 
GS and pedigree estimates have led to the conclusion 
that pedigree information may not be as useful for 
certain applications for which they have been used in 
the past (Graner et al. 1994; Plaschke et al. 1995). In 
any case, molecular-based estimates of GS will provide 
more information than is available from pedigree in- 
formation. 

Having established that molecular-based estimates 
of GS will allow plant breeders to make informed 
decisions regarding the choice of genotypes to cross, we 
must ask the question as to which assay is most appro- 
priate? Several studies have been described which ad- 
dress this question using isozymes, RFLPs and RAPDs 
(dos Santos et al. 1994; Thormann et al. 1994; Heun 
et al. 1994; Hallden et al. 1994). Heun et al. (1994) found 
that the correlation between RAPDs and isozymes 
among Avena sterilis accessions were moderately low 
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(r s = 0.36), although the overall representation of gen- 
etic relatedness was in considerable agreement. Beer 
et al. (1994) assessed genetic variation among Avena 
sterilis using morphological markers, isozymes and 
RFLPs and found a similarly low correlation 
(r 8 = 0.27). A very different situation was observed 
among Brassica species. Thormann et al. (1994) re- 
ported correlations of r 8 = 0.969 between RFLPs and 
RAPDs for a group of 18 accessions from different 
Brassica species. Dos Santos et al. (1994) also observed 
a significantly high correlation between RFLPs and 
RAPDs (r s = 0.745) using genotypes within Brassica 
oleracea, although they did observe differences between 
the RFLP and RAPD dendrograms. When Spearman 
rank correlation was used, AFLPs and RFLPs ranked 
over 70% of the pairwise comparisons in the same 
order. This may well be because both techniques are 
based on restriction site changes, the major difference is 
that PCR is used in AFLPs rather than Southern 
analysis in RFLPs. In contrast SSRs and RAPDs have 
the lowest values when compared to the other assays. 
The low correlations observed with RAPDs could be 
a reflection of the choice of primers which we have 
previously used in the construction of a linkage map 
using a population derived from two related spring 
varieties. This may well have resulted in biased esti- 
mates of GS, which in turn has affected the ranking 
order of genotypes. For example, the lowest GS 
was between two spring types (Volga and Beka; 
GS = 0.84), whereas the lowest GS estimates for the 
other assays were between two winter types (even when 
Rondo was removed from all the data sets this still 
holds true). 

The lack of correlation between SSRs and the other 
assays may not be fully unexpected, considering the 
high levels of polymorphism between pairwise com- 
parisons. Powell et al. (1996) reported that SSRs were 
well-correlated with AFLPs and RFLPs at the inter- 
species level, however at the intraspecies level the cor- 
relation disappeared, emphasising the uniqueness of 
the SSR assay. Thus, while SSR analysis appeared to be 
the most polymorphic assay system, it did not seem to 
be particularly useful for assessing genetic relationships 
among cultivars. RFLPs were particularly valuable for 
assessing genetic relationships, but required several 
probe and enzyme combinations to discriminate 
between accessions. Both RFLPs and SSRs require 
an initial investment in terms of probe or sequence 
information, and according to Vos et al. (1995) the 
ideal fingerprinting assay should require no prior se- 
quence knowledge. While only AFLPs and RAPDs 
meet these requirements, the lack of comparative in- 
formation at each assayed locus (due to dominance) 
precludes an accurate assessment of true genetic rela- 
tionships. 
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haploid [Gk. haploos, single]: Having only one set of chromosomes 
(n), in contrast to diploid (In). 

hardwood: A name commonly applied to the wood of a dicot tree. 

Hardy-Weinberg law: The mathematical expression of the relation- 
ship between the relative frequencies of two or more alleles in a pop- 
ulation; it demonstrates that the frequencies of alleles and genotypes 
will remain constant in a random-mating population in the absence 
of inbreeding, selection, or other evolutionary forces, 
haustorium, pi. haustoria [L. haustus, from haurire, to drink, draw]: 
A projection of fungal hypha that functions as a penetrating and ab- 
sorbing organ. 

heartwood: Nonliving and commonly dark-colored wood in which 
no water transport occurs; it is surrounded by sapwood. 
heliotropism [Gk. Helios, sun]: See solar tracking, 
hemicellulose (hem'i-sel'u-los): A polysaccharide resembling cellu- 
lose but more soluble and less ordered; found particularly in cell 
walls. 

herb [L. herba, grass]: A nonwoody seed plant with a relatively short- 
lived aerial portion. 

herbaceous: An adjective referring to nonwoody plants, 
herbarium: A collection of dried and pressed plant specimens, 
herbivorous: Feeding upon plants. 

heredity [L. heredis, heir]: The transmission of characteristics from 
parent to offspring through the gametes. 

hermaphrodite [Gk. for Hermes and Aphrodite]: An organism pos- 
sessing both male and female reproductive organs, 
hetero- [Gk. heteros, different]: Prefix meaning "other" or "differ- 
ent." 

heterocyst [Gk. heteros, different, + cysfzs, a bag]: A transparent, 
thick-walled, nitrogen-fixing cell that forms in the filaments of cer- 
tain cyanobacteria. 

heteroecious (hef er-e'shus) [Gk. heteros, different, + oikos, house]: 
As in some rust fungi, requiring two different host species to com- 
plete the life cycle. 

heterogamy [Gk. heteros, other, + gamos, union or reproduction]: Re- 
production involving two types of gametes. 

heterokaryotic [Gk. heteros, other, + karyon, kernel]: In fungi, having 
two or more genetically distinct types of nuclei within the same my- 
celium. 

heteromorphic [Gk. heteros, different, + morphe, form]: A term used 
to describe a life history in which the haploid and diploid generations 
are dissimilar in form. 

heterosis [Gk. heterosis, alteration]: Hybrid vigor, the superiority of 
the hybrid over either parent in any measurable character, 
heterosporous: Having two kinds of spores, designated as micro- 
spores and megaspores. 

heterothallic [Gk. heteros, different, + thallus, sprout]: A term used 
to describe a species, the haploid individuals of which are self -sterile 
or self-incompatible; two compatible strains or individuals are re- 
quired for sexual reproduction to take place. 

heterotroph [Gk. heteros, other, + trophos, feeder]: An organism that 
cannot manufacture organic compounds and so must feed on organic 
materials that have originated in other plants and animals; see also 
autotroph. 

heterozygous: Having two different alleles at the same locus on ho- 
mologous chromosomes. 

Hill reaction: The oxygen evolution and photoreduction of an artifi- 



cial electron acceptor by a chloroplast preparation in the absence of 
carbon dioxide. 

hilum [L. hilum, a trifle]: (1) Scar left on seed after separation of seed 
from funiculus; (2) the part of a starch grain around which the starch 
is laid down in more or less concentric layers, 
histone: The group of five basic proteins associated with the chromo- 
somes of all eukaryotic cells. 

holdfast: (1) Basal part of a multicellular alga that attaches it to a 
solid object; may be unicellular or composed of a mass of tissue; (2) 
cuplike structures at the tips of some tendrils, by means of which 
they become attached. 

homeo-, homo- [Gk. homos, same, similar]: Prefix meaning "similar" 
or "same." 

homeostasis (ho'me-6-sta'sis) [Gk. homos, similar, + stasis, standing]: 
The maintaining of a relatively stable internal physiological environ- 
ment within an organism, or a steady-state equilibrium in a popula- 
tion or ecosystem. Homeostasis usually involves feedback mecha- 
nisms. 

homokaryotic [Gk. homos, same, + karyon, kernel]: In fungi, having 
nuclei with the same genetic makeup within a mycelium, 
homologous chromosomes: Chromosomes that associate in pairs in 
the first stage of meiosis; each member of the pair is derived from a 
different parent. Homologous chromosomes are also called homolo- 
gues. 

homology [Gk. homologia, agreement]: A condition indicative of the 
same phylogenetic, or evolutionary, origin, but not necessarily the 
same in present structure and or function, 
homosporous: Having only one kind of spore, 
homothallic [Gk. homos, same, + thallus, sprout]: A term used to de- 
scribe a species in which the individuals are self-fertile, 
homozygous: Having identical alleles at the same locus on homolo- 
gous chromosomes. 

hormogonium, pi. hormogonia: A portion of a filament of a cyano- 
bacterium that becomes detached and grows into a new filament, 
hormone [Gk. hormaein, to excite]: A chemical substance produced 
usually in minute amounts in one part of an organism, from which it 
is transported to another part of that organism on which it has a spe- 
cific effect. 

host: An organism on or in which a parasite lives. 
* humus: Decomposing organic matter in the soU. 
hyaloplasm: See cytoplasmic ground substance, 
hybrid: Offspring of two parents that differ in one or more heritable 
characteristics; offspring of two different varieties or of two different 
species. 

hybridization: The formation of offspring between unlike parents, 
hybrid vigor: See heterosis. 

hydrocarbon [Gk. hydro, water, + L. carbo, charcoal]: An organic 
compound that consists only of hydrogen and carbon atoms, 
hydrogen bond: A weak bond between a hydrogen atom attached to 
one oxygen or nitrogen atom and another oxygen or nitrogen atom, 
hydrolysis [Gk. hydro, water, 4- lysis, loosening]: Splitting of one 
molecule into two by addition of the H + and OH" ions of water, 
hydrophyte [Gk. hydro, water, 4- phyton, a plant]: A plant that de- 
pends on an abundant supply of moisture or that ^rows wholly or 
partly submerged in water. 

hydroxyl group: An OH" group; a negatively charged ion formed by 
the dissassociation of a water molecule. 
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are no *flagella. They are *chemo- 
organotrophic, and can grow in the pres- 
ence or absence of air. They are epiphytic 
{see epiphyte) or *parasitic in plants and 
arthropods, and can cause disease (e.g. 
stubborn disease of citrus trees). 

Spirulina A genus of filamentous *cyano- 
bacteria (section III) in which the filaments 
are characteristically spiral in shape. They 
are capable of gliding motility. They are 
found in freshwater and marine habitats, 
mostly in warmer regions; and are also 
found in hot springs. Spirulina is collected 
and used as food in the region of Lake 
Chad, Africa. 

Splachnum (order *Funariales) A genus of 
mosses which are characteristically found 
growing on or associated with animal dung. 
The leaves are soft, lanceolate to obovate, 
with a thin nerve. The *capsule is held erect 
on a long *seta. Splachnum is a small genus 
found mainly in the northern hemisphere. 
S. ampullaceum and S. sphaericum occur in 
Britain. 5. ampullaceum forms dense, light- 
green tufts on cattle dung. The capsules, 
borne on red setae, are often numerous and 
are characteristic in shape, resembling mini- 
ature Greek amphorae. They are found 
mainly in the upland regions of the west 
and north of Britain, particularly on wet 
♦heaths and *moors, and in *bogs. 

spleenwort See Asplenium. 

spodic horizon A subsurface soil *horizon 
in which organic matter together with alu-. 
minium and often iron compounds have ac- ■ - 
cumulated amorphously. It is a *diagnostic 
horizon in the *USDA Soil Taxonomy. 

Spodosols An order of the *USDA Soil 
Taxonomy describing soils in which sub-. : 
surface *horizons contain amorphous mat- 
erials comprising organic matter and ' 
compounds of aluminium and often iron 
that have accumulated illuvially {see 
luviation). Such soils form in acid material, | ■ 
mainly coarse in texture, in humid cool to 
temperate climates, often beneath con-; ; |; 
iferous forests. A; 

Spondias (family *Anacardiaceae) Ajj; 
genus of trees, most of which have 'pinnate,:^ 




leaves. The fruit is a slightly angled 'drupe 
with a big round or angled stone, edible in 
several species (hog plum). There are 10 
species occurring in Indo-Malaysia and 
tropical America. 

spontaneous mutation In genetics, a natur- 
ally occurring 'mutation, as opposed to 
one artificially induced by chemicals or ir- 
radiation. Usually such mutations are due 
to errors in the normal functioning of cel- 
lular *enzymes. 

sporangiolum A 'sporangium within which 
there is only a single *spore or a small num- 
ber of spores. 

sporangiophore A specialized *hypha which 
bears a *sporangium. 

sporangiospore A ♦spore formed within a 
♦sporangium. 

sporangium A sac-like structure within 
which fungal *spores are formed; spores are 
liberated on rupture of the sporangium 
wall. 

spore 1 (mycol.) A microscopic structure 
which functions in reproduction and dis- 
persal. A spore does not contain an ♦em- 
bryo and is thus distinct from a ♦seed 
Many different types of spore are produced 
by *fungi. 2 (bacteriol.) A differentiated 
cell which may function as a 'propagule or 
as a resistant structure that allows the or- 
ganism to survive adverse environmental 
conditions. See also heterospory; homo- 
spory; megaspore; microspore; and pollen. 

spore mother cell A ♦diploid plant cell that 
gives rise to 4 *haploid ♦spores during 
*meiosis. - 

spore print The pattern of ♦spores ob- 
tained when the 'cap of a fungal ♦fruit 
body is placed, ♦gills or pores down, on a 
sheet of paper and left for a period of time. 

Sporobolomycetaceae (class ♦Blasto- 
mycetes) A family of *imperfect ♦yeasts 
characterized by the production of ♦ballis- 
tospores. These yeasts are *saprotrophic, 
and found on plant material. 

sporocarp A structure the primary func- 
tion of which is the production of ♦spores; a 
♦fruit body. The terra is applied particularly 
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to the spore-producing structure of water 
ferns, which consists of a ♦sorus completely 
enclosed by the *indusium. 

sporodochium A cushion-shaped mass of 
fungal ♦tissue densely covered with ♦conid- 
lophores. 

serogenesis The production of *spores in 
plants. In ♦mosses and 'liverworts the pro- 
cess occurs in a 'sporogonium; in ♦fungi it 
occurs in a ♦sporophore; in other spore- 
forming plants it occurs in a *sporangium. 
sporogonium In mosses (*Musci) and 
liverworts ('Hepaticae), the ♦sporophyte 
generation that develops after sexual re- 
production and produces ♦spores. 

sporogony A process of reproduction in 
which a ♦zygote undergoes multiple fission. 

sporophore A structure upon which 
'spores are borne directly. 

sporophyl] A leaf in the *Pteridophyta that 
bears 'sporangia. The term is most often 
applied to the small, scale-like leaves in the 
cones of Lycopodium and Selaginella, 
which bear ♦sporangia on their upper sides! 

sporophyte The 'spore-producing ♦dip- 
loid generation in the life cycle of plants In 
higher plants, such as *angiosperms and 
gymnosperms, the sporophyte is the dom- 
inant generation, forming the conspicuous 
plant. In lower plants, such as ♦mosses 
♦liverworts, ♦ferns, and 'fungi, the *gam- 
etophyte is the dominant and conspicuous • 
generation. See also alternation of genera- 
tions. 

sport A sudden deviation from type- a 
•mutation. 

sporuJation The process of 'spore forma- 
tion. 

spraing A disease of potatoes, which may 
be caused by either of 2 different viruses, or 
by adverse growth conditions. 'Tubers 
from infected plants show characteristic 
crescent-shaped brown marks in the flesh 
when cut. The viruses are transmitted by 
eelworms or by soil *fungi. 

spray Wood trimmings, sold in faggots for 
kindling. 
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